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Abstract 

The relativistic fluid is a highly successful model used to describe the dynamics of many- 
particle, relativistic systems. It takes as input basic physics from microscopic scales and yields 
as output predictions of bulk, macroscopic motion. By inverting the process, an understanding 
of bulk features can lead to insight into physics on the microscopic scale. Relativistic fluids 
have been used to model systems as "small" as heavy ions in collisions, and as large as the 
universe itself, with "intermediate" sized objects like neutron stars being considered along 
the way. The purpose of this review is to discuss the mathematical and theoretical physics 
underpinnings of the relativistic (multiple) fluid model. We focus on the variational principle 
approach championed by Brandon Carter and his collaborators, in which a crucial element 
is to distinguish the momenta that are conjugate to the particle number density currents. 
This approach differs from the "standard" text-book derivation of the equations of motion 
from the divergence of the stress-energy tensor, in that one explicitly obtains the relativistic 
Euler equation as an "integrability" condition on the relativistic vorticity. We discuss the 
conservation laws and the equations of motion in detail, and provide a number of (in our 
opinion) interesting and relevant applications of the general theory. 



1 Introduction 



1.1 Setting the stage 

If one does a search on the topic of relativistic fluids on any of the major physics article databases 
one is overwhelmed by the number of "hits" . This is a manifestation of the importance that 
the fluid model has long had for modern physics and engineering. For relativistic physics, in 
particular, the fluid model is essential. After all, many-particle astrophysical and cosmological 
systems are the best sources of detectable effects associated with General Relativity. Two obvious 
examples, the expansion of the universe and oscillations of neutron stars, indicate the vast range 
of scales on which relativistic fluid models are relevant. A particularly exciting context for general 
relativistic fluids today is their use in the modelling of sources of gravitational waves. This includes 
the compact binary inspiral problem, either of two neutron stars or a neutron star and a black 
hole, the collapse of stellar cores during supernovae, or various neutron star instabilities. One 
should also not forget the use of special relativistic fluids in modelling collisions of heavy nuclei, 
astrophysical jets, and gamma-ray burst sources. 

This review provides an introduction to the modelling of fluids in General Relativity. Our 
main target audience is graduate students with a need for an understanding of relativistic fluid 
dynamics. Hence, we have made an effort to keep the presentation pedagogical. The article 
will (hopefully) also be useful to researchers who work in areas outside of General Relativity and 
gravitation per se (eg. a nuclear physicist who develops neutron star equations of state), but who 
require a working knowledge of relativistic fluid dynamics. 

Throughout most of the article we will assume that General Relativity is the proper description 
of gravity. Although not too severe, this is a restriction since the problem of fluids in other 
theories of gravity has interesting aspects. As we hope that the article will be used by students 
and researchers who are not necessarily experts in General Relativity and techniques of differential 
geometry, we have included some discussion of the mathematical tools required to build models of 
relativistic objects. Even though our summary is not a proper introduction to General Relativity 
we have tried to define the tools that are required for the discussion that follows. Hopefully 
our description is sufficiently self-contained to provide a less experienced reader with a working 
understanding of (at least some of) the mathematics involved. In particular, the reader will find 
an extended discussion of the covariant and Lie derivatives. This is natural since many important 
properties of fluids, both relativistic and non-relativistic, can be established and understood by the 
use of parallel transport and Lie-dragging. But it is vital to appreciate the distinctions between 
the two. 

Ideally, the reader should have some familiarity with standard fluid dynamics, eg. at the level 
of the discussion in Landau & Lifshitz [62| . basic thermodynamics |91| . and the mathematics of 
action principles and how they are used to generate equations of motion |61| . Having stated 
this, it is clear that we are facing a facing a real challenge. We are trying to introduce a topic 
on which numerous books have been written (eg. 1021 IH7I |H1 1117] L and which requires an 
understanding of much of theoretical physics. Yet, one can argue that an article of this kind is 
timely. In particular, there has recently been exciting developments for multi-constituent systems, 
such as superfluid/superconducting neutron star cores. 1 Much of this work has been guided by the 
geometric approach to fluid dynamics championed by Carter |14M16llT%] . This provides a powerful 
framework that makes extensions to multi-fluid situations quite intuitive. A typical example of a 
phenomenon that arises naturally is the so-called entrainment effect, which plays a crucial role in 
a superfluid neutron star core. Given the potential for future applications of this formalism, we 

1 In this article we use "superfluid" to refer to any system which has the ability to flow without friction. In this 
sense, superfluids and superconductors are viewed in the same way. When we wish to distinguish charge carrying 
fluids, we will call them superconductors. 
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have opted to base much of our description on the work of Carter and his colleagues. 

Even though the subject of relativistic fluids is far from new, a number of issues remain to be 
resolved. The most obvious shortcoming of the present theory concerns dissipative effects. As we 
will discuss, dissipative effects are (at least in principle) easy to incorporate in Newtonian theory 
but the extension to General Relativity is problematic (see, for instance, Hiscock and Lindblom 
Following early work by Eckart |H6) . a significant effort was made by Israel and Stewart 
jnillM] and Carter ^]EJ- Incorporation of dissipation is still an active enterprise, and of key 
importance for future gravitational-wave asteroseismology which requires detailed estimates of the 
role of viscosity in suppressing possible instabilities. 

1.2 A brief history of fluids 

The two fluids air and water are essential to human survival. This obvious fact implies a basic 
human need to divine their innermost secrets. Homo Sapiens have always needed to anticipate air 
and water behaviour under a myriad of circumstances, such as those that concern water supply, 
weather, and travel. The essential importance of fluids for survival, and that they can be exploited 
to enhance survival, implies that the study of fluids probably reaches as far back into antiquity 
as the human race itself. Unfortunately, our historical record of this ever-ongoing study is not so 
great that we can reach very far accurately. 

A wonderful account (now in Dover print) is "A History and Philosophy of Fluid Mechanics," 
by G. A. Tokaty 105 . He points out that while early cultures may not have had universities, 
government sponsored labs, or privately funded centres pursuing fluids research (nor the Living 
Reviews archive on which to communicate results!), there was certainly some collective understand- 
ing. After all, there is a clear connection between the viability of early civilizations and their access 
to water. For example, we have the societies associated with the Yellow and Yangtze rivers in 
China, the Ganges in India, the Volga in Russia, the Thames in England, and the Seine in France, 
to name just a few. We should also not forget the Babylonians and their amazing technological 
(irrigation) achievements in the land between the Tigris and Euphrates, and the Egyptians, whose 
intimacy with the flooding of the Nile is well documented. In North America, we have the so- 
called Mississippians, who left behind their mound-building accomplishments. For example, the 
Cahokians (in Collinsville, Illinois) constructed Monk's Mound, the largest precolumbian earthen 
structure in existence that is (see http://en.wikipedia.org/wiki/Monk' s_Mound) "... over 100 
feet tall, 1000 feet long, and 800 feet wide (larger at its base than the Great Pyramid of Giza)." 

In terms of ocean and sea travel, we know that the maritime ability of the Mediterranean people 
was a main mechanism for ensuring cultural and economic growth and societal stability. The finely- 
tuned skills of the Polynesians in the South Pacific allowed them to travel great distances, perhaps 
reaching as far as South America itself, and certainly making it to the "most remote spot on the 
Earth," Easter Island. Apparently, they were remarkably adept at reading the smallest of signs 
— water color, views of weather on the horizon, subtleties of wind patterns, floating objects, birds, 
etc — as indication of nearby land masses. Finally, the harsh climate of the North Atlantic was 
overcome by the highly accomplished Nordic sailors, whose skills allowed them to reach several 
sites in North America. Perhaps it would be appropriate to think of these early explorers as adept 
geophysical fluid dynamicists/oceanographers? 

Many great scientists are associated with the study of fluids. Lost are the names of those 
individuals who, almost 400,000 years ago, carved "aerodynamically correct" wooden spears. 
Also lost are those who developed boomerangs and fin-stabilized arrows. Among those not lost is 
Archimedes, the Greek Mathematician (287-212 BC), who provided a mathematical expression for 
the buoyant force on bodies. Earlier, Thales of Miletus (624-546 BC) asked the simple question: 
What is air and water? His question is profound since it represents a clear departure from the 
main, myth-based modes of inquiry at that time. Tokaty ranks Hero of Alexandria as one of the 
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1700s FLUID DYNAMICS - 

Isaac Newton 

(1642-1727) 

Daniel Bernoulli 

(1700-1782] 

1800s VISCOSITY - 

Claude Louis Navier 

(1785-1836) 

George Gabriel Stokes 

[1819-1903] 

1900s RELATIVITY - 

Albert Einstein 

(1879-1955) 

1950s SUPERFLUIDS — 

Heike Kamerlingh-Onnes 

(1853-1926) 

Lev Davidovich Landau 

(1908-1968] 

MULTIFLUID MODELS 

Isaak Markovich Khalatnikov 

[191 9-] 

Brandon Carter 

[1942-) 



Archimedes 

(287-212 BC] 
- EARLY EMPIRICAL STUDIES 

Leonardo da Vinci 

[1425-1519) 
-VARIATIONAL METHODS 

Leonhard Euler 

[1 707-83) 

Louis de Lagrange 

[1736-1813] 



Sophus Lie 



- VORTICITY 

Hans Ertel 

[1904-1995) 

Andre Lichnerowicz 

(1915-1998) 

IRREVERSIBLE THERMODYNAMICS 

Lars Onsager 

[1903-1976] 

llya Prigogine 

(1917-2003) 

- RELATIVISTIC DISSIPATION 

Carl Henry Eckart 

(1902-1973) 

Werner Israel 

(1931) 

John M. Stewart 

(1947) 



Figure 1: A "timeline" focussed on the topics covered in this review, including chemists, engineers, 
mathematicians, philosophers, and physicists who have contributed to the development of non- 
relativistic fluids, their relativistic counterparts, multi-fluid versions of both, and exotic phenomena 
such as superfluidity. 



great, early contributors. Hero (c. 10-70) was a Greek scientist and engineer, who left behind 
many writings and drawings that, from today's perspective, indicate a good grasp of basic fluid 
mechanics. To make complete account of individual contributions to our present understanding 
of fluid dynamics is, of course, impossible. Yet it is useful to list some of the contributors to the 
field. We provide a highly subjective "timeline" in Fig.^ Our list is to a large extent focussed on 
the topics covered in this review, and includes chemists, engineers, mathematicians, philosophers, 
and physicists. It recognizes those who have contributed to the development of non-relativistic 
fluids, their relativistic counterparts, multi-fluid versions of both, and exotic phenomena such as 
superfluidity. We have provided this list with the hope that the reader can use these names as 
key words in a modern, web-based literature search whenever more information is required. 

Tokaty |1U5] discusses the human propensity for destruction when it comes to our water re- 
sources. Depletion and pollution are the main offenders. He refers to a "Battle of the Fluids" 
as a struggle between their destruction and protection. His context for this discussion was the 
Cold War. He rightly points out the failure to protect our water and air resource by the two 
predominant participants — the USA and USSR. In an ironic twist, modern study of the rela- 
tivistic properties of fluids has also a "Battle of the Fluids." A self-gravitating mass can become 
absolutely unstable and collapse to a black hole, the ultimate destruction of any form of matter. 

1.3 Notation and conventions 

Throughout the article we assume the "MTW" jZH] conventions. We also generally assume ge- 
ometrized units c = G = 1, unless specifically noted otherwise. A coordinate basis will always be 
used, with spacetime indices denoted by lowercase Greek letters that range over {0,1,2,3} (time 
being the zeroth coordinate), and purely spatial indices denoted by lowercase Latin letters that 
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range over {1, 2, 3}. Unless otherwise noted, we assume the Einstein summation convention. 
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2 Physics in a Curved Spacetime 



There is an extensive literature on special and general relativity, and the spacetime-based view 
(i.e. three space and one time dimensions that form a manifold) of the laws of physics. For 
the student at any level interested in developing a working understanding we recommend Taylor 
and Wheeler 103 for an introduction, followed by Hartle's excellent recent text 0H| designed for 
students at the undergraduate level. For the more advanced students, we suggest two of the 
classics, "MTW" [23 and Weinberg or the more contemporary book by Wald |1()8) . Finally, 

let us not forget the Living Reviews archive as a premier online source of up-to-date information! 

In terms of the experimental and/or observational support for special and general relativity, we 
recommend two articles by Will that were written for the 2005 World Year of Physics celebration 
|115l 1116] , They summarize a variety of tests that have been designed to expose breakdowns 
in both theories. (We also recommend Will's popular book Was Einstein Right? [113] and his 
technical exposition Theory and Experiment in Gravitational Physics jll4j .) To date, Einstein's 
theoretical edifice is still standing! 

For special relativity, this is not surprising, given its long list of successes: explanation of the 
Michaelson-Morley result, the prediction and subsequent discovery of anti-matter, and the standard 
model of particle physics, to name a few. Will [115] offers the observation that genetic mutations 
via cosmic rays require special relativity, since otherwise muons would decay before making it to 
the surface of the Earth. On a more somber note, we may consider the Trinity site in New Mexico, 
and the tragedies of Hiroshima and Nagasaki, as reminders of E = mc 2 . 

In support of general relativity, there are Eotvos-type experiments testing the equivalence of 
inertial and gravitational mass, detection of gravitational red-shifts of photons, the passing of 
the solar system tests, confirmation of energy loss via gravitational radiation in the Hulse- Taylor 
binary pulsar, and the expansion of the universe. Incredibly, general relativity even finds a 
practical application in the GPS system: if general relativity is neglected an error of about 15 
meters results when trying to resolve the location of an object |115) . Definitely enough to make 
driving dangerous! 

The evidence is thus overwhelming that general relativity, or at least something that passes the 
same tests, is the proper description of gravity. Given this, we assume the Einstein Equivalence 
Principle, i.e. that [TTSI HTdI ITH] : 

• test bodies fall with the same acceleration independently of their internal structure or com- 
position; 

• the outcome of any local non-gravitational experiment is independent of the velocity of the 
freely-falling reference frame in which it is performed; 

• the outcome of any local non-gravitational experiment is independent of where and when in 
the universe it is performed. 

If the Equivalence Principle holds, then gravitation must be described by a metric-based theory 
[115] . This means that 

1. spacetime is endowed with a symmetric metric, 

2. the trajectories of freely falling bodies are geodesies of that metric, and 

3. in local freely falling reference frames, the non-gravitational laws of physics are those of 
special relativity. 

For our present purposes this is very good news. The availability of a metric means that we 
can develop the theory without requiring much of the differential geometry edifice that would be 
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needed in a more general case. We will develop the description of relativistic fluids with this in 
mind. Readers that find our approach too "pedestrian" may want to consult the recent article by 
Gourgoulhon |46) . which serves as a useful complement to our description. 



2.1 The Metric and Spacetime Curvature 

Our strategy is to provide a "working understanding" of the mathematical objects that enter the 
Einstein equations of general relativity. We assume that the metric is the fundamental "field" 
of gravity. For four-dimensional spacetime it determines the distance between two spacetime 
points along a given curve, which can generally be written as a one parameter function with, 
say, components cc A '(A). As we will see, once a notion of parallel transport is established, the 
metric also encodes information about the curvature of its spacetime, which is taken to be of the 
pseudo-Riemannian type, meaning that the signature of the metric is — h H — h (cf. Eq. below). 

In a coordinate basis, which we will assume throughout this review, the metric is denoted by 
9y,v = gvy,- The symmetry implies that there are in general ten independent components (modulo 
the freedom to set arbitrarily four components that is inherited from coordinate transformations; 
cf. Eqs. (JSJl and © below). The spacetime version of the Pythagorean theorem takes the form 

ds 2 = g^dafda? , (1) 

and in a local set of Minkowski coordinates {t, x, y, z} (i.e. in a local inertial frame) it looks like 

ds 2 = - {dtf + (dxf + (dyf + {dzf . (2) 

This illustrates the — h ++ signature. The inverse metric g^ v is such that 

S^V = 8» v , (3) 

where 8^ v is the unit tensor. The metric is also used to raise and lower spacetime indices, i.e. if we 
let denote a contravariant vector, then its associated covariant vector (also known as a covector 
or one-form) is obtained as 

= g^V &V» = g^Vu . (4) 

We can now consider three different classes of curves: timelike, null, and spacelike. A vector is 
said to be timelike if g^V^V" < 0, null if g tlu V' i V v = 0, and spacelike if g^VV > 0. We can 
naturally define timelike, null, and spacelike curves in terms of the congruence of tangent vectors 
that they generate. A particularly useful timelike curve for fluids is one that is parameterized by 
the so-called proper time, i.e. x^ir) where 

dr 2 = -ds 2 . (5) 

The tangent to such a curve has unit magnitude; specifically, 

dr M 

U" = ^ , (6) 

and thus 

u v dx^dx v ds 1 

Under a coordinate transformation x^ — > x^ , contravariant vectors transform as 

V = ?V (8) 
dx v y ' 
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and covariant vectors as 

— 3t v 

v.- £,11. m 

Tensors with a greater rank (i.e. a greater number of indices), transform similarly by acting linearly 
on each index using the above two rules. 

When integrating, as we need to when we discuss conservation laws for fluids, we must be 
careful to have an appropriate measure that ensures the coordinate invariance of the integration. 
In the context of three-dimensional Euclidean space the measure is referred to as the Jacobian. 
For spacetime, we use the so-called volume form e M „ pT . It is completely antisymmetric, and for 
four-dimensional spacetime, it has only one independent component, which is 

£0123 = V^9 and e 0123 = -L= , (10) 



where g is the determinant of the metric. The minus sign is required under the square root because 
of the metric signature. By contrast, for three-dimensional Euclidean space (i.e. when considering 
the fluid equations in the Newtonian limit) we have 

£123 = V9 and e 123 = ±- , (11) 

where g is the determinant of the three-dimensional space metric. A general identity that is 
extremely useful is |108| 

^... NN+1 ...^^ jVj+ ^ Mn = { _ 1 y {n _ m5 [^ v ^. 5 ^] vn (12) 

where s is the number of minus signs in the metric (eg. s = 1 here). As an application consider 
three-dimensional, Euclidean space. Then s = and n = 3. Using lower-case Latin indices we 
see (for j = 1 in Eq. (15) ) 

e mij e mkl = 8 l k 8h - V k 5\ . (13) 

The general identity will be used in our discussion of the variational principle, while the three- 
dimensional version is often needed in discussions of Newtonian vorticity. 

2.2 Parallel Transport and the Covariant Derivative 

In order to have a generally covariant prescription for fluids, i.e. written in terms of spacetime 
tensors, we must have a notion of derivative that is itself covariant. For example, when 
acts on a vector V v a rank- two tensor of mixed indices must result: 

„ Qx a dx v 

- fe^ 1 " ■ < 14 > 

The ordinary partial derivative does not work because under a general coordinate transformation 

dV V dx a dx v dVP dx a d 2 x u 



dx M dx^ dxP dx a dx» dx a dxP 



V p . (15) 



The second term spoils the general covariance, since it vanishes only for the restricted set of 
rectilinear transformations 

= a»x v + bP , (16) 

where a% and 6 M are constants. 

For both physical and mathematical reasons, one expects a covariant derivative to be defined in 
terms of a limit. This is, however, a bit problematic. In three-dimensional Euclidean space limits 
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Figure 2: A schematic illustration of two possible versions of parallel transport. In the first case 
(a) a vector is transported along great circles on the sphere locally maintaining the same angle with 
the path. If the contour is closed, the final orientation of the vector will differ from the original 
one. In case (b) the sphere is considered to be embedded in a three-dimensional Euclidean space, 
and the vector on the sphere results from projection. In this case, the vector returns to the original 
orientation for a closed contour. 

can be defined uniquely in that vectors can be moved around without their lengths and directions 
changing, for instance, via the use of Cartesian coordinates, the {i, j, k} set of basis vectors, and the 
usual dot product. Given these limits those corresponding to more general curvilinear coordinates 
can be established. The same is not true of curved spaces and/or spacetimes because they do not 
have an a priori notion of parallel transport. 

Consider the classic example of a vector on the surface of a sphere (illustrated in Fig. |2J) . Take 
that vector and move it along some great circle from the equator to the North pole in such a way 
as to always keep the vector pointing along the circle. Pick a different great circle, and without 
allowing the vector to rotate, by forcing it to maintain the same angle with the locally straight 
portion of the great circle that it happens to be on, move it back to the equator. Finally, move 
the vector in a similar way along the equator until it gets back to its starting point. The vector's 
spatial orientation will be different from its original direction, and the difference is directly related 
to the particular path that the vector followed. 

On the other hand, we could consider that the sphere is embedded in a three-dimensional 
Euclidean space, and that the two-dimensional vector on the sphere results from projection of a 
three-dimensional vector. Move the projection so that its higher-dimensional counterpart always 
maintains the same orientation with respect to its original direction in the embedding space. When 
the projection returns to its starting place it will have exactly the same orientation as it started 
with, see Fig. [5] It is thus clear that a derivative operation that depends on comparing a vector 
at one point to that of a nearby point is not unique, because it depends on the choice of parallel 
transport. 

Pauli 83 notes that Levi-Civita is the first to have formulated the concept of parallel 
"displacement," with Weyl [112] generalizing it to manifolds that do not have a metric. The point 
of view expounded in the books of Weyl and Pauli is that parallel transport is best defined as a 
mapping of the "totality of all vectors" that "originate" at one point of a manifold with the totality 
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at another point. (In modern texts, one will find discussions based on fiber bundles.) Pauli points 
out that we cannot simply require equality of vector components as the mapping. 

Let us examine the parallel transport of the force-free, point particle velocity in Euclidean 
three-dimensional space as a means for motivating the form of the mapping. As the velocity is 
constant, we know that the curve traced out by the particle will be a straight line. In fact, we 
can turn this around and say that the velocity parallel transports itself because the path traced 
out is a geodesic (i.e. the straightest possible curve allowed by Euclidean space). In our analysis 
we will borrow liberally from the excellent discussion of Lovelock and Rund [7T|. Their text is 
comprehensive yet readable for one who is not well-versed with differential geometry. Finally, 
wc note that this analysis will be of use later when we obtain the Newtonian limit of the general 
relativistic equations, but in an arbitrary coordinate basis. 

We are all well aware that the points on the curve traced out by the particle are describable, in 
Cartesian coordinates, by three functions x l (t) where t is the standard Newtonian time. Likewise, 
we know the tangent vector at each point of the curve is given by the velocity components v l (t) = 
dx 1, /dt, and that the force-free condition is equivalent to 

a\t) = ^ = =S> v\t) = const . (17) 

Hence, the velocity components v l (0) at the point x l (0) are equal to those at any other point along 
the curve, say v l (T) at x l (T), and so we could simply take v z (0) = v l (T) as the mapping. But as 
Pauli warns, we only need to reconsider this example using spherical coordinates to see that the 
velocity components {r, 9, <fr} must necessarily change as they undergo parallel transport along a 
straight-line path (assuming the particle does not pass through the origin) . The question is what 
should be used in place of component equality? The answer follows once we find a curvilinear 
coordinate version of dv z /dt = 0. 

What we need is a new "time" derivative D/di, that yields a generally covariant statement 

where the u*(i) = dx 1 /dt are the velocity components in a curvilinear system of coordinates. Let 
us consider an infinitesimal displacement of the velocity from x l (t) to x l (t + St). In Cartesian 
coordinates we know that 

Sv* = ^-6x3 = . (19) 

dxi K ' 

Consider now a coordinate transformation to the curvilinear coordinate system x l , the inverse 
being x % = x i (x J ). Given that 

v% = ^ ( 2 °) 

■^—v 3 ) v k . (21) 
dt ydxi dx k dx k dx 3 J V ; 

The metric c/,. for our curvilinear coordinate system is obtainable from 



we can write 

dv l ( dx 1 dv J 



where 



_ dx k dx 1 

9ij = -^-[Tzz;9ki , (22) 
ox ox J 



'■>•> = \ I • (23) 
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Differentiating Eq. (12211 with respect to x, and permuting indices, we can show that 
d 2 x h dx l 1 



dx^x 3 dx 2 

where 



k9hi = -x (fljfcj + g jk ,i ~ ffij,fc) = 3u { /* } , (24) 



With some tensor aigebra we find 



*W - |^ ■ ( 25 ) 



dw* _ 9x l Dw' 

dt" ~ ft? ~dT ' ( 6) 



where _ 

DtJ 1 ■ / dv 



One can verify that the operator D/di is covariant with respect to general transformations of 
curvilinear coordinates. 

We now identify our generally covariant derivative (dropping the overline) as 

Similarly, the covariant derivative of a covector is 

dv ■ 

^'•< {■)'■>. r, :i . (29) 

One extends the covariant derivative to higher rank tensors by adding to the partial derivative 
each term that results by acting linearly on each index with { * } using the two rules given above. 

By relying on our understanding of the force-free point particle, we have built a notion of 
parallel transport that is consistent with our intuition based on equality of components in Cartesian 
coordinates. We can now expand this intuition and see how the vector components in a curvilinear 
coordinate system must change under an infinitesimal, parallel displacement from x l (t) to x l (t+St). 
Setting Eq. (|2"T|) to zero, and noting that v t 5t = Sx l , implies 

5v i =^Sx* = -{ h ' j }v k 6xl . (30) 

In general relativity we assume that under an infinitesimal parallel transport from a spacetime 
point x M (A) on a given curve to a nearby point ^(A + 5X) on the same curve, the components of 
a vector will change in an analogous way, namely 

SV? = ^-5x u = -Vt v V p 5x u , (31) 
II q x v py \ ' 

where 

5x^ = -tt-SX . (32) 
dA 

Weyl (112] refers to the symbol as the "components of the affine relationship," but we will use 
the modern terminology and call it the connection. In the language of Weyl and Pauli, this is the 
mapping that we were looking for. 

For Euclidean space, we can verify that the metric satisfies 

Vigjk = (33) 
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for a general, curvilinear coordinate system. The metric is thus said to be "compatible" with the 
covariant derivative. We impose metric compatibility in general relativity. This results in the 
so-called Christoffel symbol for the connection, defined as 

T U = 7^.9 P<T (SW + AW - 9»v,<t) ■ (34) 

The rules for the covariant derivative of a contravariant vector and a covector are the same as in 
Eqs. I|28[) and l|29[) . except that all indices are for full spacetime. 



2.3 The Lie Derivative and Spacetime Symmetries 

Another important tool for measuring changes in tensors from point to point in spacetime is the 
Lie derivative. It requires a vector field, but no connection, and is a more natural definition in 
the sense that it does not even require a metric. It yields a tensor of the same type and rank 
as the tensor on which the derivative operated (unlike the covariant derivative, which increases 
the rank by one). The Lie derivative is as important for Newtonian, non-relativistic fluids as for 
relativistic ones (a fact which needs to be continually emphasized as it has not yet permeated the 
fluid literature for chemists, engineers, and physicists). For instance, the classic papers on the 
Chandrasekhar-Friedman-Schutz instability ^2 LH] m rotating stars are great illustrations of the 
use of the Lie derivative in Newtonian physics. We recommend the book by Schutz [SS] for a 
complete discussion and derivation of the Lie derivative and its role in Newtonian fluid dynamics 
(see also the recent series of papers by Carter and Chamel ^OBOH^)- We will adapt here the 
coordinate-based discussion of Schouten |51|. as it may be more readily understood by those not 
well- versed in differential geometry. 

In a first course on classical mechanics, when students encounter rotations, they are introduced 
to the idea of active and passive transformations. An active transformation would be to fix the 
origin and axis-orientations of a given coordinate system with respect to some external observer, 
and then move an object from one point to another point of the same coordinate system. A passive 
transformation would be to place an object so that it remains fixed with respect to some external 
observer, and then induce a rotation of the object with respect to a given coordinate system, by 
rotating the coordinate system itself with respect to the external observer. We will derive the Lie 
derivative of a vector by first performing an active transformation and then following this with 
a passive transformation and finding how the final vector differs from its original form. In the 
language of differential geometry, we will first "push-forward" the vector, and then subject it to a 
"pull-back" . 

In the active, push-forward sense we imagine that there are two spacetime points connected by 
a smooth curve x M (A). Let the first point be at A = 0, and the second, nearby point at A = e, 
i.e. £ M (e); that is, 

a^=a^(e) » a# + e , (35) 

where Xq = x^(0) and 

e = ^- (36) 

dA A=0 

is the tangent to the curve at A = 0. In the passive, pull-back sense we imagine that the coordinate 
system itself is changed to x^ = x fJ, (x"), but in the very special form 



= x" - e ^ ■ (37) 

It is in this last step that the Lie derivative differs from the covariant derivative. In fact, if we 
insert Eq. ^ into gZJ we find the result x 1 ^ = Xq . This is called "Lie-dragging" of the coordinate 
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Figure 3: A schematic illustration of the Lie derivative. The coordinate system is dragged along 
with the flow, and one can imagine an observer "taking derivatives" as he/she moves with the flow 
(see the discussion in the text). 



frame, meaning that the coordinates at A = are carried along so that at A = e, and in the new 
coordinate system, the coordinate labels take the same numerical values. 

As an interesting aside it is worth noting that Arnold UJ, only a little whimsically, refers to this 
construction as the "fisherman's derivative" . He imagines a fisherman sitting in a boat on a river, 
"taking derivatives" as the boat moves with the current. Playing on this imagery, the covariant 
derivative is cast with the high-test Zebco parallel transport fishing pole, the Lie derivative with 
the Shimano, Lie-dragging ultra-light. Let us now see how Lie-dragging reels in vectors. 

For some given vector field that takes values V*(A), say, along the curve, we write 



= V>" (0) 



(38) 



for the value of at A = and 



Vf = 7"(e) 



(39) 



for the value at A = e. Because the two points Xq and x^ are infinitesimally close e << 1, and we 
can thus write 

dV* 1 



(40) 



A=0 



for the value of at the nearby point and in the same coordinate system. However, in the new 
coordinate system at the nearby point, we find 



dx v 



(41) 



A=0 



The Lie derivative now is defined to be 



VP - yp 
LcV* = lim —i 

£U _____ yV __ 

~ dx v dx v 
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(42) 



where we have dropped the "0" subscript and the last equality follows easily by noting = T^. 

We introduced in Eq. the effect of parallel transport on vector components. By contrast, 
the Lie-dragging of a vector causes its components to change as 

6Vg = C^V e . (43) 

We see that if C^V^ = 0, then the components of the vector do not change as the vector is Lie- 
dragged. Suppose now that represents a vector field and that there exists a corresponding 
congruence of curves with tangent given by £ M . If the components of the vector field do not change 
under Lie-dragging we can show that this implies a symmetry, meaning that a coordinate system 
can be found such that the vector components will not depend on one of the coordinates. This is 
a potentially very powerful statement. 

Let £^ represent the tangent to the curves drawn out by, say, the fj, = a coordinate. Then we 
can write x a (X) = A which means 

e = ■ (44) 
If the Lie derivative of V 11 with respect to vanishes we find 

r^r— = V v ^- = . (45) 
s dx v dx w v ; 

Using this in Eq. (|40|) implies Vf = Vq, that is to say, the vector field V fi (x v ) does not depend 
on the x a coordinate. Generally speaking, every £ M that exists that causes the Lie derivative of a 
vector to vanish represents a symmetry. 



2.4 Spacetime Curvature 

The main message of the previous two subsections is that one must have an a priori idea of 
how vectors and higher rank tensors are moved from point to point in spacetime. Another 
manifestation of the complexity associated with carrying tensors about in spacetime is that the 
covariant derivative does not commute. For a vector we find 



VpV^V - V CT V P V^ = R^upaV , (46) 

where R^ vpa is the Riemann tensor. It is obtained from 

7->m — _ pM i pM p-r _ pM pr (a>j\ 

vpa — J- va .p J- vp<a i T p L vo 1 tct 1 vp ■ V ti I 

Closely associated are the Ricci tensor R va = R av and scalar R that are defined by the contractions 

Rva = R^ufia , R = Rva ■ (48) 

We will later also have need for the Einstein tensor, which is given by 

= R^ - ]^Rgp U ■ (49) 

It is such that V pG^ u vanishes identically (this is known as the Bianchi Identity). 

A more intuitive understanding of the Riemann tensor is obtained by seeing how its presence 
leads to a path-dependence in the changes that a vector experiences as it moves from point to 
point in spacetime. Such a situation is known as a "non-integrability" condition, because the 
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result depends on the whole path and not just the initial and final points. That is, it is unlike a 
total derivative which can be integrated and thus depends on only the lower and upper limits of 
the integration. Geometrically we say that the spacetime is curved, which is why the Riemann 
tensor is also known as the curvature tensor. 

To illustrate the meaning of the curvature tensor, let us suppose that we are given a surface 
that can be parameterized by the two parameters A and r\. Points that live on this surface will 
have coordinate labels x^{X 1 rj). We want to consider an infinitesimally small "parallelogram" 
whose four corners (moving counterclockwise with the first corner at the lower left) are given by 
x tJ, (X,rj), x^{X,n + 8rj), x tl (X + SX,n + Sn), and x t± (X + 5X, rj). Generally speaking, any "movement" 
towards the right of the parallelogram is effected by varying 77, and that towards the top results 
by varying A. The plan is to take a vector V M (A, 77) at the lower-left corner x fl (X,r]), parallel 
transport it along a A = constant curve to the lower-right corner at x^(X,n + $77) where it will 
have the components V fl (X, n + 6rf), and end up by parallel transporting V^ 1 at x M (A, 77 + Sn) along 
an ry = constant curve to the upper-right corner at x tl (X + <5A, n + Sn). We will call this path I 
and denote the final component values of the vector as . We repeat the same process except 
that the path will go from the lower-left to the upper-left and then on to the upper-right corner. 
We will call this path II and denote the final component values as . 

Recalling Eq. I|31[) as the definition of parallel transport, we first of all have 

V»(X, n + Sn)^ V»(X, 77) + 5,V((X, rj) = V(X, rj) ~ T^VS^r (50) 

and 

V»(X + 5X, n) « W (A, 77) + 5 X V[(X, n) = n) - T^ p V v 5 x x" , (51) 

where 

5 n x» sa x" ( A, 77 + Srj) - x» (A, rj) , 8 x x^ « of* (A + 8X, n) - x^ ( A, 77) . (52) 
Next, we need 

« ViX^ + S^+dxV^X^ + Sr]) (53) 

and 

V£ « V»(X + 6X,r 1 ) + 6 v V l [ l (X + 5X,T 1 ) . (54) 
Working things out, we find that the difference between the two paths is 

AV = Vf - V& = R» 1/pa V5 x x°5 ri x p , (55) 
which follows because SxS^x^ 1 — 6 v 5xx^, i.e. we have closed the parallelogram. 
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3 The Stress-Energy-Momentum Tensor and the Einstein 
Equations 

Any discussion of relativistic physics must include the stress-energy-momentum tensor T^. It is 
as important for general relativity as G^ v in that it enters the Einstein equations in as direct a 
way as possible, i.e. 

Gfj, v = 8irT^ . (56) 

Misner, Thorne, and Wheeler [7S] refer to T^ v as "...a machine that contains a knowledge of the 
energy density, momentum density, and stress as measured by any and all observers at that event." 

Without an a priori, physically-based specification for T^ u , solutions to the Einstein equations 
are devoid of physical content, a point which has been emphasized, for instance, by Geroch and 
Horowitz (in 49 ). Unfortunately, the following algorithm for producing "solutions" has been 
much abused: (i) specify the form of the metric, typically by imposing some type of symmetry, 
or symmetries, (ii) work out the components of G^ v based on this metric, (iii) define the energy 
density to be Goo and the pressure to be Gn, say, and thereby "solve" those two equations, and 
(iv) based on the "solutions" for the energy density and pressure solve the remaining Einstein 
equations. The problem is that this algorithm is little more than a mathematical game. It is 
only by sheer luck that it will generate a physically viable solution for a non-vacuum spacetime. 
As such, the strategy is antithetical to the raison d'etre of gravitational- wave astrophysics, which 
is to use gravitational-wave data as a probe of all the wonderful microphysics, say, in the cores 
of neutron stars. Much effort is currently going into taking given microphysics and combining it 
with the Einstein equations to model gravitational-wave emission from realistic neutron stars. To 
achieve this aim, we need an appreciation of the stress-energy tensor and how it is obtained from 
microphysics. 

Those who are familiar with Newtonian fluids will be aware of the roles that total internal 
energy, particle flux, and the stress tensor play in the fluid equations. In special relativity we 
learn that in order to have spacetime covariant theories (e.g. well-behaved with respect to the 
Lorentz transformations) energy and momentum must be combined into a spacetime vector, whose 
zeroth component is the energy and the spatial components give the momentum. The fluid stress 
must also be incorporated into a spacetime object, hence the necessity for T^ v . Because the 
Einstein tensor's covariant divergence vanishes identically, we must have also V^T^^ = (which 
we will later see happens automatically once the fluid field equations are satisfied). 

To understand what the various components of T M „ mean physically we will write them in terms 
of projections into the timelike and spacelike directions associated with a given observer. In order 
to project a tensor index along the observer's timelike direction we contract that index with the 
observer's unit four- velocity U^. A projection of an index into spacelike directions perpendicular 
to the timelike direction defined by (see (21] for the idea from a "3+1" point of view, or ^2] 
from the "brane" point of view) is realized via the operator _L^, defined as 

J_£= + EW„ , EW P = -1. (57) 

Any tensor index that has been "hit" with the projection operator will be perpendicular to the 
timelike direction associated with t/ M in the sense that U v = 0. Fig. 0Jis a local view of both 
projections of a vector for an observer with unit four- velocity U^. More general tensors are 
projected by acting with or _L(^ on each index separately (i.e. multi- linearly) . 

The energy density p as perceived by the observer is (see Eckart j^E] for one of the earliest 
discussions) 

p = T^U^U" , (58) 
= - J_g U v T pv (59) 
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J W*' = timelike worldline 



Figure 4: The projections at point P of a vector onto the worldline defined by ?7 M and into the 
perpendicular hypersurface (obtained from the action of -L£). 

is the spatial momentum density, and the spatial stress iS M „ is 

<V =JLJJ-g T CTP . (60) 

The manifestly spatial component Sij is understood to be the i th -component of the force across 
a unit area that is perpendicular to the j th -direction. With respect to the observer, the stress- 
energy-momentum tensor can be written in full generality as the decomposition 

T„ v = P U^U U + 2U {li V u) + <V - (61) 

where 2U^V u) = U^V V + U V V^. Because U^V^ = 0, we see that the trace T = T% is 

T = S - p , (62) 

where S — S^^. We should point out that use of p for the energy density is not universal. Many 
authors prefer to use the symbol e and reserve p for the mass-density. We will later (in Sec. I6.2JI 
use the above decomposition as motivation for the simplest perfect fluid model. 
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4 Why are fluids useful models? 



The Mcrriam- Webster online dictionary (http : //www .m-w . com/ ) defines a fluid as ". . .a substance 
(as a liquid or gas) tending to flow or conform to the outline of its container" when taken as a noun 
and ". . .having particles that easily move and change their relative position without a separation of 
the mass and that easily yield to pressure: capable of flowing" when taken as an adjective. The best 
model of physics is the Standard Model which is ultimately the description of the "substance" that 
will make up our fluids. The substance of the Standard Model consists of remarkably few classes of 
elementary particles: leptons, quarks, and so-called "force" carriers (gauge- vector bosons). Each 
elementary particle is quantum mechanical, but the Einstein equations require explicit trajectories. 
Moreover, cosmology and neutron stars are basically many particle systems and, even forgetting 
about quantum mechanics, it is not practical to track each and every "particle" that makes them 
up. Regardless of whether these are elementary (leptons, quarks, etc) or collections of elementary 
particles (eg. stars in galaxies and galaxies in cosmology). The fluid model is such that the inherent 
quantum mechanical behaviour, and the existence of many particles are averaged over in such a 
way that it can be implemented consistently in the Einstein equations. 

Central to the model is the notion of "fluid particle," also known as a fluid element. It is 
an imaginary, local "box" that is infinitesimal with respect to the system en masse and yet large 
enough to contain a large number of particles (eg. an Avogadro's number of particles). This is 
illustrated in Fig. EI In order for the fluid model to work we require M >> N >> 1 and D » L. 
Strictly speaking, our model has L infinitesimal, M — > oo, but with the total number of particles 
remaining finite. The explicit trajectories that enter the Einstein equations are those of the fluid 
elements, not the much smaller (generally fundamental) particles that are "confined," on average, 
to the elements. Hence, when we speak later of the fluid velocity, we mean the velocity of fluid 
elements. 




Figure 5: An object with a characteristic size D is modelled as a fluid that contains M fluid 
elements. From inside the object we magnify a generic fluid element of characteristic size L. In 
order for the fluid model to work we require M » N » 1 and D » L. 

In this sense, the use of the phrase "fluid particle" is very apt. For instance, each fluid clement 
will trace out a timelike trajectory in spacetime. This is illustrated in Fig.0for a number of fluid 
elements. An object like a neutron star is a collection of worldlines that fill out continuously a 
portion of spacetime. In fact, we will see later that the relativistic Euler equation is little more 
than an "integrability" condition that guarantees that this filling (or fibration) of spacetime can be 
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performed. The dual picture to this is to consider the family of three-dimensional hypersurfaccs 
that are pierced by the worldlines at given instants of time, as illustrated in Fig.0 The integrability 
condition in this case will guarantee that the family of hypersurfaces continuously fill out a portion 
of spacetime. In this view, a fluid is a so-called three-brane (see |18| for a general discussion of 
branes). In fact the method used in Sec. [SI to derive the relativistic fluid equations is based on 
thinking of a fluid as living in a three-dimensional "matter" space (i.e. the left-hand-side of Fig. EJ. 

Once one understands how to build a fluid model using the matter space, it is straight-forward 
to extend the technique to single fluids with several constituents, as in Sec. and multiple fluid 
systems, as in Sec. EH An example of the former would be a fluid with one species of particles 
at a non-zero temperature, i.e. non-zero entropy, that does not allow for heat conduction relative 
to the particles. (Of course, entropy does flow through spacetime.) The latter example can 
be obtained by relaxing the constraint of no heat conduction. In this case the particles and 
the entropy are both considered to be fluids that are dynamically independent, meaning that the 
entropy will have a four-velocity that is generally different from that of the particles. There is 
thus an associated collection of fluid elements for the particles and another for the entropy. At 
each point of spacetime that the system occupies there will be two fluid elements, in other words, 
there are two matter spaces (cf. Sec. 110(1 . Perhaps the most important consequence of this is that 
there can be a relative flow of the entropy with respect to the particles (which in general leads to 
the so-called entrainmcnt effect). The canonical example of such a two fluid model is superfluid 
He 4 [Sj]. 
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5 A Primer on Thermodynamics and Equations of State 



Fluids consist of many fluid elements, and each fluid element consists of many particles. The 
state of matter in a given fluid element is determined thermodynamically |91|. meaning that only 
a few parameters are tracked as the fluid element evolves. Generally, not all the thermodynamic 
variables arc independent, being connected through the so-called equation of state. The number 
of independent variables can also be reduced if the system has an overall additivity property. As 
this is a very instructive example, we will now illustrate such additivity in detail. 



5.1 Fundamental, or Euler, Relation 

Consider the standard form of the combined first and second laws for a simple, single-species 
system: 

dE = TdS - pdV + pdN . (63) 
This follows because there is an equation of state, meaning that E = E(S, V, N) where 



T=^ 
dS 



dE 

V,N dV 



_ dE_ 
S ,N ' ^~9N 



(64) 

s,v 



The total energy E, entropy S, volume V, and particle number N are said to be extensive if when 
S, V, and N are doubled, say, then E will also double. Conversely, the temperature T, pressure 
p, and chemical potential p are called intensive if they do not change their values when V, N, and 
S are doubled. This is the additivity property and we will now show that it implies an Euler 
relation (also known as the "fundamental relation" |91p among the thermodynamic variables. 

Let a tilde represent the change in thermodynamic variables when S, V, and N are all increased 
by the same amount A, i.e. 

S = XS , V = XV , N = XN . (65) 
Taking E to be extensive then means 

E(S, V, N) = XE(S, V, N) . (66) 

Of course we have for the intensive variables 

f = T, p = p , fi, = ii . (67) 

Now, 

dE = XdE + EdX 

= TdS - pdV + fidN 

= X (TdS - pdV + ndN) + (TS -pV + pN) dX , (68) 

and therefore we find the Euler relation 

E = TS — pV + /j,N . (69) 

If we let p = E/V denote the total energy density, s = S/V the total entropy density, and n — N/V 
the total particle number density, then 

p + p = Ts + [in . (70) 

The nicest feature of an extensive system is that the number of parameters required for a 
complete specification of the thermodynamic state can be reduced by one, and in such a way that 
only intensive thermodynamic variables remain. To see this, let A = l/V, in which case 

S = s , V = 1 , N = n . (71) 
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The re-scaled energy 
p = p(s, n) since 



The first law thus becomes 



becomes just the total energy density, i.e. E = E/V = p, and moreover 
p = E(S, V, N) = E{S/V, 1, N/V) = E(s, n) . (72) 



oi- 



dS = TdS - pdV + pdN = Tds + pdn 



dp = Tds + pdn 



This implies 



T = 

ds 



_ dp 
dn 



The Euler relation Eq. (|7()JI then yields the pressure as 



P 



dp 

ds 



dp 
dn 



(72) 

(73) 
(74) 
(75) 

(76) 



We can think of a given relation p(s,n) as the equation of state, to be determined in the flat, 
tangent space at each point of the manifold, or, physically, a small enough region across which the 
changes in the gravitational field are negligible, but also large enough to contain a large number of 
particles. For example, for a neutron star Glendenning has reasoned that the relative change 
in the metric over the size of a nucleon with respect to the change over the entire star is about 
10 -19 , and thus one must consider many internucleon spacings before a substantial change in the 
metric occurs. In other words, it is sufficient to determine the properties of matter in special 
relativity, neglecting effects due to spacetime curvature. The equation of state is the major link 
between the microphysics that governs the local fluid behaviour and global quantities (such as the 
mass and radius of a star). 

In what follows we will use a thermodynamic formulation that satisfies the fundamental scaling 
relation, meaning that the local thermodynamic state (modulo entrainment) is a function of the 
variables N/V, S/V, etcetera. This is in contrast to the fluid formulation of "MTW" 75 1. In 
their approach one fixes from the outset the total number of particles N, meaning that one simply 
sets dN = in the first law of thermodynamics. Thus without imposing any scaling relation, one 
can write 

dp = d (E/V) = Tds+-(p + p- Ts) dn . (77) 
n 

This is consistent with our starting point for fluids, because we assume that the extensive variables 
associated with a fluid element do not change as the fluid element moves through spacetime. 
However, we feel that the use of scaling is necessary in that the fully conservative, or non-dissipative, 
fluid formalism presented below can be adapted to non-conservative, or dissipative, situations where 
diV = cannot be imposed. 



5.2 Prom Microscopic Models to Fluid Equation of State 

Let us now briefly discuss how an equation of state is constructed. For simplicity, we focus on a one- 
parameter system, with that parameter being the particle number density. The equation of state 
will then be of the form p = p(n). In many-body physics (such as studied in condensed matter, 
nuclear, and particle physics) one can in principle construct the quantum mechanical particle 
number density ^qm, stress-energy- momentum tensor T^ 1 ^, and associated conserved particle 
number density current «q M (starting with some fundamental Lagrangian, say; cf. |ll)9M45lllll)) ). 
But unlike in quantum field theory in curved spacetime ^2Jj one assumes that the matter exists 
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in an infinite Minkowski spacetime (cf. the discussion following Eq. JSJ)). If the reader likes, the 
application of T^ M at a spacetime point means that T^ 1 has been determined with respect to a 
flat tangent space at that point. 

Once T^ M is obtained, and after (quantum mechanical and statistical) expectation values with 
respect to the system's (quantum and statistical) states are taken, one defines the energy density 
as 

p = u^(T^), (78) 

where 

At sufficiently small temperatures, p will just be a function of the number density of particles n at 
the spacetime point in question, i.e. p — p(n). Similarly, the pressure is obtained as 

P=l((T QM %)+p) (80) 

and it will also be a function of n. 

One must be very careful to distinguish T^ M from T^ v . The former describes the states of 
elementary particles with respect to a fluid element, whereas the latter describes the states of fluid 
elements with respect to the system. Comer and Joynt |.'i2j have shown how this line of reasoning 
applies to the two-fluid case. 
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6 An Overview of the Perfect Fluid 



There are many different ways of constructing general relativistic fluid equations. Our purpose 
here is not to review all possible methods, but rather to focus on a couple: (i) an "off-the-shelve" 
consistency analysis for the simplest fluid a la Eckart to establish some key ideas, and then (ii) a 
more powerful method based on an action principle that varies fluid element world lines. The ideas 
behind this variational approach can be traced back to Taub |1()2) (see also |HHJ)- Our description 
of the method relies heavily on the work of Brandon Carter, his students, and collaborators [1611331 
1341 1251 l2lfl 1631 1551 186) . We prefer this approach as it utilizes as much as possible the tools of the 
trade of relativistic fields, i.e. no special tricks or devices will be required (unlike even in the case of 
our "off-the-shelve" approach). One's footing is then always made sure by well-grounded, action- 
based derivations. As Carter has always made clear: when there are multiple fluids, of both the 
charged and uncharged variety, it is essential to distinguish the fluid momenta from the velocities. 
In particular, in order to make the geometrical and physical content of the equations transparent. 
A well-posed action is, of course, perfect for systematically constructing the momenta. 

6.1 Rates-of-Change and Eulerian Versus Lagrangian Observers 

The key geometric difference between generally covariant Newtonian fluids and their general rel- 
ativistic counterparts is that the former have an a priori notion of time 1201 l2*T] . Newtonian 
fluids also have an a priori notion of space (which can be seen explicitly in the Newtonian covariant 
derivative introduced earlier; cf. the discussion in Such a structure has clear advantages for 

evolution problems, where one needs to be unambiguous about the rate-of-change of the system. 
But once a problem requires, say, electromagnetism, then the a priori Newtonian time is at odds 
with the full spacetime covariance of electromagnetic fields. Fortunately, for spacetime covariant 
theories there is the so-called "3+1" formalism (see, for instance, [OH]) that allows one to define 
"rates-of-change" in an unambiguous manner, by introducing a family of spacelike hypersurfaces 
(the "3") given as the level surfaces of a spacetime scalar (the "1"). 

Something that Newtonian and relativistic fluids have in common is that there are preferred 
frames for measuring changes — those that are attached to the fluid elements. In the parlance 
of hydrodynamics, one refers to Lagrangian and Eulerian frames, or observers. A Newtonian 
Eulerian observer is one who sits at a fixed point in space, and watches fluid elements pass by, 
all the while taking measurements of their densities, velocities, etc. at the given location. In 
contrast, a Lagrangian observer rides along with a particular fluid element and records changes of 
that element as it moves through space and time. A relativistic Lagrangian observer is the same, 
but the relativistic Eulerian observer is more complicated to define. Smarr and York |22J define 
such an observer as one who would follow along a worldline that remains everywhere orthogonal 
to the family of spacclikc hypersurfaces. 

The existence of a preferred frame for a one fluid system can be used to great advantage. In 
the very next sub-section we will use an "off-the-shelve" analysis that exploits the fact of one 
preferred frame to derive the standard perfect fluid equations. Later, we will use Eulerian and 
Lagrangian variations to build an action principle for the single and multiple fluid systems. These 
same variations can also be used as the foundation for a linearized perturbation analysis of neutron 
stars 59 . As we will see, the use of Lagrangian variations is absolutely essential for establishing 
instabilities in rotating fluids |41II42) . Finally, we point out that multiple fluid systems can have 
as many notions of Lagrangian observer as there are fluids in the system. 
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6.2 The Single, Perfect Fluid Problem: "Off-the-shelve" Consistency 
Analysis 

We earlier took the components of a general stress-energy-momentum tensor and projected them 
onto the axes of a coordinate system carried by an observer moving with four- velocity U^. As 
mentioned above, the simplest fluid is one for which there is only one four- velocity . Hence, 
there is a preferred frame defined by and if we want the observer to sit in this frame we can 
simply take If 1 = u^. With respect to the fluid there will be no momentum flux, i.e. = 0. Since 
we use a fully spacetime covariant formulation, i.e. there are only spacetime indices, the resulting 
stress-energy-momentum tensor will transform properly under general coordinate transformations, 
and hence can be used for any observer. 

The spatial stress is a two-index, symmetric tensor, and the only objects that can be used to 
carry the indices are the four- velocity and the metric g^ v . Furthermore, because the spatial 
stress must also be symmetric, the only possibility is a linear combination of g^ v and u^u". Given 
that u^S^v — 0, we find 

S^v = + UpUv) . (81) 

If we identify the pressure as p — 5/3 then 

Tp V = (p + p) u^Uy + pg^y , (82) 

which is the well-established result for a perfect fluid. 

Given a relation p = p(p), there are four independent fluid variables. Because of this the 
equations of motion are often understood to be V^T'V = 0, which follows immediately from the 
Einstein equations and the fact that V^G^ 1 ,, = 0. To simplify matters, we take as equation of state 
a relation of the form p = p(n) where n is the particle number density. The chemical potential is 
then given by 

dp = — — An = pdn , (83) 
an 

and we see from the Euler relation Eq. (|70|l that 

pn=p + p . (84) 



Let us now get rid of the free index of V^T^ 1 ,, = 
and second, by projecting it with (letting — 
have the identity 

v m ko = o =► 

Contracting with u v and using this identity gives 



in two ways: first, by contracting it with u v 
u M ). Recalling the fact that = -1 we 

u v V^u y = . (85) 



tt"V M p + (p+p)V M u'* = ■ ( 86 ) 
The definition of the chemical potential p and the Euler relation allow us to rewrite this as 

pu^V^n + pnV^ => V M ?i M = , (87) 

where n M = nu^ . Projection of the free index using _L£ leads to 

D»P = -{p + pW . (88) 

where =-L£ V p is a purely spatial derivative and a M = u u V „u M is the acceleration. This is 
reminiscent of the familiar Euler equation for Newtonian fluids. 
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/„ - 2n^ v , (93) 



However, we should not be too quick to think that this is the only way to understand V ^T^v = 
0. There is an alternative form that makes the perfect fluid have much in common with vacuum 
electromagnetism. If we define 

Mm = /-"V ' ( 89 ) 
and note that u M du M = (because u^u^ — — 1), then 

dp = -// M dn^ . (90) 

The stress-energy-momentum tensor can now be written in the form 

T" v = pS" v + n^u u . (91) 

We have here our first encounter with the momentum fi^ that is conjugate to the particle number 
density current n M . Its importance will become clearer as this review develops, particularly when 
we discuss the two-fluid case. If we now project onto the free index of V^T^ 1 ,, = using _L^, as 
before, we find 

fu + (V„n") n v = , (92) 

where the force density /„ is 
and the vorticity is defined as 

c*V = v [m/M = \ - V^) . (94) 

Contracting Eq. I|92|l with n v implies (since w M „ = — a-v^) that 

= (95) 

and as a consequence 

/„ = 2n**w^ = . (96) 

The vorticity two-form lu^ u has emerged quite naturally as an essential ingredient of the fluid 
dynamics |67l 1161 1101 I57| . Those who are familiar with Newtonian fluids should be inspired 
by this, as the vorticity is often used to establish theorems on fluid behaviour (for instance the 
Kelvin-Hclmholtz Theorem and is at the heart of turbulence modelling |88| . 

To demonstrate the role of lo^ u as the vorticity, consider a small region of the fluid where the 
time direction t^, in local Cartesian coordinates, is adjusted to be the same as that of the fluid 
four-velocity so that = t^ = (1,0,0,0). Eq. (|96|) and the antisymmetry then imply that ui^ 
can only have purely spatial components. Because the rank of u>^ v is two, there are two "nulling" 
vectors, meaning their contraction with either index of u; M „ yields zero (a condition which is true 
also of vacuum electromagnetism). We have arranged already that t 1 * be one such vector. By 
a suitable rotation of the coordinate system the other one can be taken as = (0,0,0, 1), thus 
implying that the only non-zero component of uj^ is u) xy . "MTW" 17^1 points out that such a 
two-form can be pictured geometrically as a collection of oriented worldtubes, whose walls lie in the 
x = constant and y = constant planes. Any contraction of a vector with a two-form that does not 
yield zero implies that the vector pierces the walls of the worldtubes. But when the contraction is 
zero, as in Eq. (|96|l . the vector does not pierce the walls. This is illustrated in Fig.® where the red 
circles indicate the orientation of each world-tube. The individual fluid element four-velocities lie 
in the centers of the world-tubes. Finally, consider the dashed, closed contour in Fig. If that 
contour is attached to fluid-element worldlines, then the number of worldtubes contained within 
the contour will not change because the worldlines cannot pierce the walls of the worldtubes. This 
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Figure 6: A local, geometrical view of the Euler equation as an integrability condition of the 
vorticity for a single-constituent perfect fluid. 



is essentially the Kelvin- Helmholtz Theorem on conservation of vorticity. From this we learn that 
the Euler equation is an integrability condition which ensures that the vorticity two-surfaces mesh 
together to fill out spacetime. 

As we have just seen, the form Eq. (|96|) of the equations of motion can be used to discuss the 
conservation of vorticity in an elegant way. It can also be used as the basis for a derivation of 
other known theorems in fluid mechanics. To illustrate this, let us derive a generalized form of 
Bernoulli's theorem. Let us assume that the flow is invariant with respect to transport by some 
vector field k p . That is, we have 



(V p /i ff - V CT At p )fc ff = V P (£;V CT ) 



(97) 



Here one may consider two particular situations. If k a is taken to be the four-velocity, then the 
scalar k a ' \i a represents the "energy per particle" . If instead k a represents an axial generator of 
rotation, then the scalar will correspond to an angular momentum. For the purposes of the present 
discussion we can leave k" unspecified, but it is still useful to keep these possibilities in mind. Now 
contract the equation of motion Eq. (|§fj|) with k a and assume that the conservation law Eq. 
holds. Then it is easy to show that we have 



n p V p (ii a k a ) = V p {n p n„k°) = 



(98) 



In other words, we have shown that n p /i cr /c <T is a conserved quantity. 

Given that we have just inferred the equations of motion from the identity that V^T^V = 0, we 
now emphatically state that while the equations are correct the reasoning is severely limited. In 
fact, from a field theory point of view it is completely wrong! The proper way to think about the 
identity is that the equations of motion are satisfied first, which then guarantees that V^T^ 1 ,, = 0. 
There is no clearer way to understand this than to study the multi-fluid case: Then the vanishing 
of the covariant divergence represents only four equations, whereas the multi-fluid problem clearly 
requires more information (as there are more velocities that must be determined). We have 
reached the end of the road as far as the "off-the-shelf" strategy is concerned, and now move on 
to an action-based derivation of the fluid equations of motion. 
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7 Setting the Context: The Point Particle 



The simplest physics problem, i.e. the point particle, has always served as a guide to deep principles 
that are used in much harder problems. We have used it already to motivate parallel transport 
as the foundation for the covariant derivative. Let us call upon the point particle again to set the 
context for the action-based derivation of the fluid field equations. We will simplify the discussion 
by considering only motion in one dimension. We assure the reader that we have good reasons, 
and ask for patience while we remind him/her of what may be very basic facts. 
Early on we learn that an action appropriate for the point particle is 

I = [ 1 dt T = / ' dt (\mx 2 \ , (99) 



where m is the mass and T the kinetic energy. A first-order variation of the action with respect 
to x{t) yields 



51 = -J dt (mx)Sx + {mx5x)\\l . (100) 

If this is all the physics to be incorporated, i.e. if there are no forces acting on the particle, then 
we impose d'Alembert's principle of least action [SJ, which states that those trajectories x[t) that 
make the action stationary, i.e. SI = 0, are those that yield the true motion. We see from the 
above that functions x(t) that satisfy the boundary conditions 

Sx(ti) = = 5x(tf) , (101) 

and the equation of motion 

mx = , (102) 

will indeed make 81 = 0. This same reasoning applies in the substantially more difficult fluid 
actions that will be considered later. 

But, of course, forces need to be included. First on the list are the so-called conservative forces, 
describable by a potential V(x), which are placed into the action according to: 

/ = jf 1 dtL(x, x) = jf ' dt Qmi 2 - V{x) \ , (103) 

where L = T — V is the so-called Lagrangian. The variation now leads to 

81 = -J 1 dt (mx+ 5x+ (m±5x)\ll . (104) 

Assuming no externally applied forces, d'Alembert's principle yields the equation of motion 

dV 

rax + — = . (105) 
ox 

An alternative way to write this is to introduce the momentum p (not to be confused with the 
fluid pressure introduced earlier) defined as 

dL . . 

P = ttt = fnx , (106) 
ox 

in which case 

dV 

P + - = 0. (107) 
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In the most honest applications, one has the obligation to incorporate dissipative, i.e. non- 
conservative, forces. Unfortunately, dissipative forces Fd cannot be put into action principles. 
Fortunately, Newton's second law is of great guidance, since it states 



mx + — — = F d , 



(108) 



when conservative and dissipative forces act. A crucial observation of Eq. (|108(l is that the 
"kinetic" (mx = p) and conservative (dV/dx) forces, which enter the left-hand side, do follow from 
the action, i.e. 



When there are no dissipative forces acting, the action principle gives us the appropriate equation 
of motion. When there are dissipative forces, the action defines for us the kinetic and conservative 
force terms that are to be balanced by the dissipative terms. It also defines for us the momentum. 

We should emphasize that this way of using the action to define the kinetic and conservative 
pieces of the equation of motion, as well as the momentum, can also be used in a context where the 
system experiences an externally applied force F ext . The force can be conservative or dissipative, 
and will enter the equation of motion in the same way as Fd did above. That is 



Like a dissipative force, the main effect of the external force can be to siphon kinetic energy from 
the system. Of course, whether a force is considered to be external or not depends on the a priori 
definition of the system. To summarize: The variational argument leads to equations of motion 
of the form 




(109) 



T" — Fd + F ext ■ 



(110) 




SI .. dV . dV 



= Conservative 

7^ Dissipative and/or External Forces 



(111) 



dx dx 
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8 The "Pull-back" Formalism for a Single Fluid 



In this section the equations of motion and the stress-energy- momentum tensor for a one-component, 
general relativistic fluid are obtained from an action principle. Specifically a so-called "pull-back" 
approach (see, for instance, 133 021 0^) is used to construct a Lagrangian displacement of the 
number density four-current n^, whose magnitude n is the particle number density. This will 
form the basis for the variations of the fundamental fluid variables in the action principle. 

As there is only one species of particle considered here, n M is conserved, meaning that once a 
number of particles N is assigned to a particular fluid element, then that number is the same at 
each point of the fluid element's worldline. This would correspond to attaching a given number 
of particles (i.e. Ni, N2, etc.) to each of the worldlines in Fig. [7] Mathematically, one can write 
this as a standard particle-flux conservation equation: 

= . (112) 

For reasons that will become clear shortly, it is useful to rewrite this conservation law in what may 
(if taken out context) seem a rather obscure way. We introduce the dual n v \ T to n M , i.e. 

n v xr = e vX r^ , n» = ^ vXT n v \ T , (113) 

such that 

n 2 = ^n„ XT n» XT , (114) 

In Fig. we have seen that a two-form gives worldtubes. A three-form is the next higher-ranked 
object and it can be thought of, in an analogous way, as leading to boxes [75]. The key step here 
is to realize that the conservation rule is equivalent to having the three-form n v \ T be closed. In 
differential geometry this means that 

V [M n„ AT] = . (115) 

This can be shown to be equal to Eq. I|112fl by contracting with ^ vXr . 

The main reason for introducing the dual is that it is straightforward to construct a particle 
number density three-form that is automatically closed, since the conservation of the particle 
number density current should not — speaking from a strict field theory point of view — be a part of 
the equations of motion, but rather should be automatically satisfied when evaluated on a solution 
of the "true" equations. 

This can be made to happen by introducing a three-dimensional "matter" space — the left-hand 
part of Fig. — which can be labelled by coordinates X A = X A (x fJ -), where A, B, C etc = 1, 2, 3. 
For each time slice in spacetime, there will be a corresponding configuration in the matter space; 
that is, as time goes forward, the weaving of worldlines in spacetime will be matched by the fluid 
particle positions in the matter space. In this sense we are "pulling-back" from spacetime to the 
matter space (cf. the previous discussion of the Lie derivative) . The three-form can be pulled-back 
to its three-dimensional matter space by using the mappings X A , as indicated in Fig. This 
allows us to construct a three-form that is automatically closed on spacetime. Hence, we let 

n vXr = N ABC (VuX a ) (V x X b ) \7 t X c , (116) 

where Nabc is completely antisymmetric in its indices and is a function only of the X A . The X A 
will also be comoving in the sense that 

n^^X A = le^ XT N BC D (V„A S ) (V A A C ) (V T X D ) \7 ^X A = . (117) 
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Figure 7: The pull-back from "fluid-particle" worldlincs in spacctime, on the right, to "fluid- 
particle" points in the three-dimensional matter space labelled by the coordinates {X 1 , X 2 , X 3 }. 
Here, the pull-back of the "/ th " (/ — 1,2, ...,n) worldline into the abstract space (for the t = 
configuration of "fluid-particle" points) is Xf = X A (0, x\) where x\ is the spatial position of the 
intersection of the worldline with the t = time slice. 

The time part of the spacetime dependence of the X A is thus somewhat ad hoc; that is, if we would 
take the flow of time to be the proper time of the worldlines, then the X A would not change, and 
hence there would be no "motion" of the fluid particles in the matter space. 

Because the matter space indices are three-dimensional and the closure condition involves four 
spacetime indices, and also the X A are scalars on spacetime (and thus two covariant differentiations 
commute), the pull-back construction does indeed produce a closed three- form: 

V [m 7V At] = V [M (N ABC V y X A V x X B V T] X c ) = . (118) 

In terms of the scalar fields X A , we now have particle number density currents that are automati- 
cally conserved. Thus, another way of viewing the pull-back construction is that the fundamental 
fluid field variables are the X A (as evaluated on spacetime). In fact, the variations of n v \ T can 
now be taken with respect to variations of the X A . 

Let us introduce the Lagrangian displacement on spacetime for the particles, to be denoted £ M . 
This is related to the variation SX A via a "push-forward" construction (which takes a variation in 
the matter space and pushes it forward to spacetime): 

SX A = - (V„X A ) . (119) 

Using the fact that 

vjx A = -v v ([v^x A ]e) 

= - (V^X A ) - (V^ V X A ) (120) 

means 

Sn^r = - (fV ff n„Ar + n CTA rV^ CT + ?wV A r + n„A CT V T n - -£tn u xr , (121) 
where C is the Lie derivative along £ M (cf. Eq. I|42|l ). We can thus infer that 

6n» - n a W a e - £ a V CT n" - n» + ^g ap Sg ap ) . (122) 
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By introducing the decomposition 



= nu^ 1 , = -1 , (123) 

we can furthermore show that 



Sn = -V CT - n ( u v u°V a t v + - [g°r + n^u?] 8g ap ) , (124) 



1 

2 

and 

5u» = (6£ + u»u p ) (u a \7 a e - eVaU p ) + ^u»u° U PSg ap . (125) 
The Lagrangian variation resulting from the Lagrangian displacement is given by 

A = S + C^ , (126) 

from which it follows that 

An^ T = , (127) 
which is entirely consistent with the pull-back construction. We also find that 



1 

2 

1 
2 

and 



Au^ = -vPu' T u p Ag ap , (128) 



Ae uX ra = -€ v xto9 w Ag^p , (129) 
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An = --(g°r + u a uP)Ag a p . (130) 

These formulae (and their Newtonian analogues) have been adroitly used by Friedman and Schutz 
in their establishment of the so-called Chandrasekhar-Friedman-Schutz (CFS) instability |27l 1411 
132]. see Sec.P 

With a general variation of the conserved four-current in hand, we can now use an action 
principle to derive the equations of motion and the stress-energy-momentum tensor. The central 
quantity in the analysis is the so-called "master" function A, which is a function of the scalar n 2 = 
—n^nP. For this single fluid system, it is such that —A corresponds to the local thermodynamic 
energy density. In the action principle, the master function is the Lagrangian density for the fluid: 

i4 



Ifiuid = / V^Ad 4 x . (131) 

An unconstrained variation of A(n 2 ) is with respect to n M and the metric gp U , and allows the 
four components of to be varied independently. It takes the form 

SA = iipSn^ + inV*^ . (132) 

where 

^ = Bn p , B = -20 . (133) 

The use of the letter B is to remind us that this is a bulk fluid effect, which is present regardless 
of the number of fluids and constituents. The momentum covector \x p is what we encountered 
earlier. It is dynamically, and thermodynamically, conjugate to n^, and its magnitude can be 
seen to be the chemical potential of the particles by recalling that A = —p. Correctly identifying 
the momentum is essential for a concise and transparent formulation of fluids, e.g. for constructing 
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proofs about vorticity, and also for coupling in dissipative terms in the equations of motion (cf. the 
discussion of Sec.EJ. For later convenience, we also introduce the momentum form defined as 

= e ^X PfJ/p ( 134 ) 

If the variation of the four-current was left unconstrained, the equations of motion for the fluid 
deduced from the above variation of A would require, incorrectly, that the momentum covector /j,^ 
should vanish in all cases. This reflects the fact that the variation of the conserved four-current 
must be constrained, meaning that not all components of n M can be treated as independent. In 
terms of the constrained Lagrangian displacement of Eq. i|122[l , a first-order variation of the master 
function A yields 

9M V 

(135) 

where the total divergence does not contribute to the field equations or stress-energy-momentum 
tensor (the divergence theorem implies that it becomes a boundary term in the action), the force 
density f v is as defined in Eq. (|9"3T) . and the generalized pressure \& is defined to be 

* = A - nV P • (136) 

We have utilized the well-known formula 75 

SV^9 = \y/=9<r&9v» ■ (137) 

At this point we can return to the view that n M is the fundamental field for the fluid. For 
vanishing metric variations, the principle of least action implies that the coefficient of as well 
as the total divergence (which yields the boundary term), must vanish. Thus, the equations of 
motion consist of the original conservation condition from Eq. (|112l) and the Euler equation 

/„ = 2n"c<v = . (138) 

We see that (j,"* 7 n\ T M £ M , as evaluated on the boundary, is fixed during the variation. 

When 6g^ v ^ 0, it is well-established in the relativity literature that the stress-energy-momentum 
tensor follows from a variation of the Lagrangian with respect to the metric, that is 

T^5 9lMV = -=<5 (V=#C) - + n^x) g Xu 5 9fW . (139) 

When the complete set of field equations is satisfied then we see from Eq. (|96l) . which applies here, 
that it is automatically true that V^T^ 1 ,, = 0. 

Let us now recall the discussion of the point particle. There we pointed out that only the 
fully conservative form of Newton's Second Law follows from an action, i.e. external or dissipative 
forces are excluded. However, we argued that a well-established form of Newton's Second Law 
is known that allows for external and/or dissipative forces, cf. Eq. I|lllfl . There is thus much 
purpose in using the particular symbol f v in Eq. (|138|) . We may take the f v to be the relativistic 
analogue of the left-hand-side of Eq. (| 1 1 1 1) in every sense. In particular, when dissipation and/or 
external "forces" act in a general relativistic setting, they are incorporated via the right-hand-side 
of Eq. ifTSg)) . 



^A) 



v 1 r 



g^ Xr n XT ,e 
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9 The Two-Constituent, Single Fluid 



Generally speaking, the total energy density p can be a function of independent parameters other 
than the particle number density n n , such as the entropy density n s , assuming that the system 
scales in the manner discussed in Sec. [5] so that only densities need enter the equation of state. 
(For later convenience we will introduce the constituent indices x, y, etc., which range over the 
two constituents {n, s}, and do not satisfy any kind of summation convention.) If there is no 
heat conduction, then this is still a single fluid problem, meaning that there is still just one unit 
flow velocity |33) . This is what we mean by a two-constituent, single fluid. We assume 
that the particle number and entropy are both conserved along the flow, in the same sense as in 
Eq. I|117|l . Associated which each parameter there is thus a conserved current density four- vector, 
i.e. = rinU^ for the particle number density and = rigM^ for the entropy density. Note that 
the ratio x s = n s /n n is comoving in the sense that 



(x a ) = . (140) 



In terms of constituent indices x, y, the associated first law can be written in the form 

dp = J2 V K dn x = - ti d < . ( 141 ) 

x— {n.s} x— {n.s} 



since p = p(n^,n1), where 
and 



n£ = n x u^ , = -g^n^ , (142) 

^=v^X , BX ^ 2 ^- ( 143 ) 



Because of Eq. (|74^l p s =T where T is the temperature. 

Given that we only have one four-velocity, the system will still just have one fluid element per 
spacetime point. But unlike before, there will be an additional conserved number, N s , that can 
be attached to each worldline, like the particle number N n of Fig. [7| In order to describe the 
worldlines we can use the same three scalars X A (x tl ) as before. But how do we get a construction 
that allows for the additional conserved number? Recall that the intersections of the worldlines 
with some hypersurface, say t = 0, is uniquely specified by the three X A (0,x l ) scalars. Each 
worldline will have also the conserved numbers N n and iV s assigned to them. Thus, the values 
of these numbers can be expressed as functions of the X A (0, x l ). But most importantly, the fact 
that each N x is conserved, means that this function specification must hold for all of spacetime, 
so that in particular the ratio x s is of the form x s (x fJ ') = x s (X A (x^)). Consequently, we now have 
a construction whereby this ratio identically satisfies Eq. I|14()J) . and the action principle remains 
a variational problem just in terms of the three X A scalars. 

The variation of the action follows like before, except now a constituent index x must be 
attached to the particle number density current and three-form: 

Kxr = eu\r„< ■ (144) 

Once again it is convenient to introduce the momentum form, now defined as 

fi£ uX = e^ A V* . (145) 

Since the X A are the same for each n* Xr , the above discussion indicates that the pull-back con- 
struction now is to be based on 

Kxr = N\ BC (V„X A ) (V A X B ) \7 T X C , (146) 
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where N\ BC is completely antisymmetric and a function of the X A . After some due deliberation, 
the reader should be convinced that the only thing required here in addition to the previous 
subsection is to attach a x index to n^\ T , n^, and n in Eqs. i|121fl . i|122fl . and l|124|l . respectively. 
If we now define the master function as 

A = -p (147) 

and the generalized pressure '5 to be 

* = a- ]T mX = a + E ( 148 ) 

x— {n,s} x— {n,s} 

then the first-order variation for A is 

5(V=^A) = l -J—g^cr + [*-A]u^)5g^-^—g\ ]T f*\?+ 



ix={n,s} 



(149) 



where 
and 

The final equations of motion are 
and 



H = 2n>* i/ , (150) 

- V [X] ■ ( 151 ) 

E /? = ' ( 152 ) 

x— {n,s} 



V„< = , (153) 
whereas the stress-energy-momentum tensor takes the familiar form 

T iw = + - AJu'V . (154) 
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10 The "Pull-Back" Formalism for Two Fluids 

Having discussed the single fluid model, and how one accounts for stratification, it is time to 
move on to the problem of modelling multi-fluid systems. We will experience for the first time 
novel effects due to the presence of a relative flow between two, interpenetrating fluids, and the 
fact that there is no longer a single, preferred rest-frame. This kind of formalism is necessary, 
for example, for the simplest model of a neutron star, since it is generally accepted that the 
inner crust is permeated by an independent neutron superfluid, and the outer core is thought to 
contain superfluid neutrons, superconducting protons, and a highly degenerate gas of electrons. 
Still unknown is the number of independent fluids that would be required for neutron stars that 
have colour- flavour-locked quark matter in the deep core [Q. The model can also be used to 
describe superfluid Helium and heat-conducting fluids, which is of importance for incorporation 
of dissipation (see Sec. 114(1 . We will focus on this example here, ft should be noted that, even 
though the particular system we concentrate on consists of only two fluids it illustrates all new 
features of a general multi-fluid system. Conceptually, the greatest step is to go from one to two 
fluids. A generalization to a system with more degrees of freedoms is straightforward. 

In keeping with the previous section, we will rely on use of the constituent index, which for all 
remaining formulas of this section will range over x, y, z = n, s. In the example that we consider 
the two fluids represent the particles (n) and the entropy (s). Once again, the number density 
four-currents, to be denoted n£, are taken to be separately conserved, meaning 

V„n£ = . (155) 

As before, we use the dual formulation, i.e. introduce the three-forms 



1 

3! 



x v\t — CiAt/^x i n x — Tn e ^ XTn uXT ■ (156) 



Also like before, the conservation rules are equivalent to the individual three-forms being closed, 

V[X At] = . (157) 

However, we need a formulation whereby such conservation obtains automatically, at least in 
principle. 

We make this happen by introducing the three-dimensional matter space, the difference being 
that we now need two such spaces. These will be labelled by coordinates X£, and we recall that 
A,B,C etc = 1,2,3. This is depicted in Fig. EI which indicates the important facts that (i) a 
given spatial point can be intersected by each fluid's worldline and (ii) the individual worldlines 
are not necessarily parallel at the intersection, i.e. the independent fluids are interpenetrating and 
can exhibit a relative flow with respect to each other. Although we have not indicated this in 
Fig. [5] (in order to keep the figure as uncluttered as possible) attached to each worldline of a given 
constituent will be a fixed number of particles Nf, N$, etc. (cf. Fig. 0. For the same reason, we 
have also not labelled (as in Fig. |7|) the "pull-backs" (represented by the arrows) from spacetime 
to the matter spaces. 

By "pulling-back" each constituent's three- form onto its respective matter space we can once 
again construct three-forms that are automatically closed on spacetime, i.e. let 

nlxr = NI BC V V X^V X X^V T X^ , (158) 

where N^ BC is completely antisymmetric in its indices and is a function of the X^. Using the 
same reasoning as in the single fluid case, the pull-back construction produces three-forms that 
are automatically closed, i.e. they satisfy Eq. I|157fl identically If we let the scalar fields X^ (as 
functions on spacetime) be the fundamental variables, they yield a representation for each particle 
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Figure 8: The pull-back from an x -constituent "fluid-particle" worldlinc in spacetime, on the 
right, to "fluid-particle" points in that x th -constituent's three-dimensional matter space labelled 
by the coordinates {X^,X^,X^}, where the constituent index x = n, s. There exist as many 
matter spaces as there are dynamically independent fluids, which for this case means two. 



number density current that is automatically conserved. The variations of the three-forms can 
now be derived by varying them with respect to the X£- 

Lagrangian displacements on spacetime for each fluid, to be denoted will now be introduced. 
They are related to the variations 8X^ via the "push-forwards" 

SXf = - (V M A^) & . (159) 

It should be clear that the analogues of Eqs. I|121|l . l|122[l . (|124|) . and (|125() for this two-fluid case are 
given by the same formulas except that each displacement and four-current will now be associated 
with a constituent index, using the decomposition 

<=n x < , u*< = -l, (160) 

where u* — g^u^. 

Associated with each constituent's Lagrangian displacement is its own Lagrangian variation. 
These are naturally defined to be 

A x = <S + £ 5x , (161) 

so that it follows that 

A X <Ar - , (162) 

as expected for the pull-back construction. Likewise, two-fluid analogues of Eqs. (11281) (|130|) exist 
which take the same form except that the constituent index is attached. However, in contrast 
to the ordinary fluid case, there are many more options to consider. For instance, we could also 
look at the Lagrangian variation of the first constituent with respect to the second constituent's 
flow, i.e. A s n n , or the other way around, i.e. A n n s . The Newtonian analogues of these Lagrangian 
displacements were essential to an analysis of instabilities in rotating superfluid neutron stars . 
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We are now in a position to construct an action principle that will yield the equations of motion 
and the stress-energy tensor. Again, the central quantity is the "master" function A, which is now 
a function of all the different scalars that can be formed from the n£, i.e. the scalars n x together 
with 

n xy = n yx = ~9^K n y ■ ( 16 3) 

In the limit where all the currents are parallel, i.e. the fluids are comoving, —A corresponds again 
to the local thermodynamic energy density. In the action principle, A is the Lagrangian density for 
the fluids. It should be noted that our choice to use only the fluid currents to form scalars implies 
that the system is "locally isotropic" in the sense that there are no a priori preferred directions, 
i.e. the fluids are equally free to move in any direction. Structures like the crust believed to exist 
near the surface of a neutron star generally would be locally anisotropic, eg. sound waves in the 
lattice move in the preferred directions associated with the lattice. 

An unconstrained variation of A with respect to the independent vectors n^. and the metric g^ 
takes the form 

SA= ]T /i*<K + ^ A!/ ( E <Ma)%,, (164) 

x— {n,s} \x— {n,s} / 

where 

M£ = 9^ (#X + A**r%) , (165) 

A xy = A y, = _9A_ ^ _ (lg6) 

The momentum covectors /z* are each dynamically, and thermodynamically, conjugate to their 
respective number density currents n£, and their magnitudes are the chemical potentials. We see 
that A* y represents the fact that each fluid momentum /z* may, in general, be given by a linear 
combination of the individual currents n£. That is, the current and momentum for a particular 
fluid do not have to be parallel. This is known as the entrainment effect. We have chosen to 
represent it by the letter A for historical reasons. When Carter first developed his formalism he 
opted for this notation, referring to the "anomaly" of having misaligned currents and momenta. 
It has since been realized that the entrainment is a key feature of most multi-fluid systems and it 
would, in fact, be anomalous to leave it out! 

In the general case, the momentum of one constituent carries along some mass current of the 
other constituents. The entrainment only vanishes in the special case where A is independent 
of n xy (x 7^ y) because then we obviously have A xy = 0. Entrainment is an observable effect 
in laboratory superfluids |89l 111)4) (eg. via flow modifications in superfluid 4 He and mixtures of 
superfluid 3 He and 4 He). In the case of neutron stars, entrainment is an essential ingredient of 
the current best explanation for the so-called glitches U2JI221- Carter has also argued that 
these "anomalous" terms are necessary for causally well-behaved heat conduction in relativistic 
fluids, and by extension necessary for building well-behaved relativistic equations that incorporate 
dissipation. 

In terms of the constrained Lagrangian displacements, a variation of A now yields 



-<?A) 



\ x— {n,s} / x— {n,s} 



(167) 
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where /* is as defined in Eq. H150JI except that the individual velocities are no longer parallel. The 
generalized pressure $ is now 

¥ = A- n >l- ( 168 ) 

x— {n,s} 

At this point we will return to the view that and n£ are the fundamental variables for the 
fluids. Because the are independent variations, the equations of motion consist of the two 
original conservation conditions of Eq. (|155|l plus two Euler type equations 

11 = . (169) 

We also find that the stress-energy tensor is 

T"„ = + • (170) 

x— {n,s} 

When the complete set of field equations is satisfied then it is automatically true that V^T^ 1 ,, = 0. 
One can also verify that is symmetric. The momentum form entering the boundary 

term is the natural extension of Eq. I|145fl to this two-fluid case. 

It must be noted that Eq. l|169|l is significantly different from the multi-constituent version of 
Eq. (|152l) . This is true even if one is solving for a static and spherically symmetric configura- 
tion, where the fluid four-velocities would all necessarily be parallel. Simply put, Eq. (|169|l still 
represents two independent equations. If one takes entropy as an independent fluid, then the 
static and spherically symmetric solutions will exhibit thermal equilibrium |35) . This explains, 
for instance, why one must specify an extra condition (eg. convective stability |lllp to solve for a 
double-constituent star with only one four-velocity. 
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11 Speeds of Sound 



Crucial to the understanding of black holes is the effect of spacetime curvature on the light-cone 
structures, that is, the totality of null vectors that emanate from each spacetime point. Crucial 
to the propagation of massless fields is the light-cone structure. In the case of fluids, it is both the 
speed of light and the speed (and/or speeds) of sound that dictate how waves propagate through 
the matter, and thus how the matter itself propagates. We give here a simple analysis that uses 
plane-wave propagation to derive the speed(s) of sound in the single-fluid, two-constituent single- 
fluid, and two-fluid cases, assuming that the metric variations vanish (see for a more rigorous 
derivation) . The analysis is local, assuming that the speed of sound is a locally defined quantity, 
and performed using local, Minkowski coordinates x M . The purpose of the analysis is to illuminate 
how the presence of various constituents and multi-fluids impact the local dynamics. 

In a small region we will assume that the configuration of the matter with no waves present 
is locally isotropic, homogeneous, and static. Thus for the background n£ = {n x , 0, 0, 0} and the 
vorticity uf^ v vanishes. For all cases the general form of the variation of the force density /* for 
each constituent is 

Sfl = 2<3 [ „<5 M * ] • (171) 
Similarly, the conservation of the n£ gives 

<9 M< S< = . (172) 

We are here using the point of view that the n£ are the fundamental fluid fields, and thus plane 
wave propagation means 

6n£ = A»e ik » x " , (173) 
where the amplitudes A£ and the wave vector fc M are constant. Eqs. i|172fl and 1173fl imply simply 

k^Sn^ = ; (174) 

i.e. the waves are "transverse" in the spacetime sense. We will give general forms for <5/i* in what 
follows, and only afterwards assume static, homogeneous, and isotropic backgrounds. 



11.1 Single Fluid Case 

Suppose that there is only one constituent, with index x = n. The master function A then depends 
only on n^. The variation in the chemical potential due to a small disturbance Sn£ is 

6^ = BIMl , (175) 

where 

B% = B n g,„ - 2—g^g vp nlni . (176) 

The equation of motion is 6f° = 0. It is not difficult to show, by using the condition of transverse 
wave propagation Eq. (1 17411 and contracting with the spatial part of the wave vector k % = g lJ kj, 
that the equation of motion reduces to 



The speed of sound is thus 



B n + B n o^2-J^< = 0. (177) 
^ 2 B dlnB n , . 

c ^¥tr 1 + 8^;- (178) 

Given that a well-constructed fluid model should have C| < 1, we see that the second term must 
be negative. This would ensure that the model is causal. 
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11.2 Two-constituent, Single Fluid Case 



Now consider the case when there are the two constituents with densities n n and n s , two conserved 
density currents n£ and n^ 1 , two chemical potential covectors /x£ and /i^, but still only one four- 
velocity u^. The master function A depends on both and meaning that 

8fi=BZ v 6fZ + X$8r%, (179) 

where 

dB x 

X 7» = - 2 Qrf9^9» P <n p y (180) 

and B^ is as in Eq. I|176fl except that each n is replaced with s. We have chosen to use the letter 
X to represent what is a true multi-constituent effect, which arises due to composition gradients 
in the system. An alternative would have been to use C since the effect is due to the presence 
of different constituents. However, in his papers Carter tends to use B s = C, referring to the 
bulk entropy contribution as "caloric" . Our chosen notation is intended to avoid confusion. It is 
also the case that the presence of the composition term has not been emphasized in previous 
work. This may be unfortunate since a composition variation is known to affect the dynamics of 
a system, eg. by giving rise to the g-modes in a neutron star |93| . 

The fact that n£ is parallel to n£ implies that it is only the magnitude of the entropy density 
current that is independent. One can show that the condition of transverse propagation, as applied 
to both currents, implies 

<5< = x B 5n£ ■ (181) 

Now, we proceed as in the previous example, noting that the equation of motion is 5f™ + Sf^ = 0, 
which reduces to 

[B n + x 2 s B s ] i - [B n cl + x 2 B s c 2 - 2x s X™} ) M< - , (182) 



where 



The speed of sound is thus 



^ 1 + §e£- <»> 



K?n 2 I 2bs 2 _ o yns 

° s ~ WTxW* ■ (184) 



11.3 Two Fluid Case 

The two-fluid problem is qualitatively different from the previous two cases, since there are now two 
independent density currents. This fact impacts the analysis in two crucial ways: (i) the master 
function A depends on n^, nf, and n„ s = n 2 n (i.e. entrainment is present), and (ii) the equations 
of motion, after taking into account the transverse flow condition, are doubled to 5/° = = 6f*. 
As we will see, the key point is that there are now two sound speeds that must be determined. 
A variation of the chemical potential covectors that leaves the metric fixed takes the form 

*A*Z - + A™) SK + + A%) Sn; , (185) 

where the complicated terms 

A™ = -g,.g vp [<n° + «] + ^nXj , (186) 
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/ <9S X dB y dA xy \ 

= A^g^ - g^g vp [q^T« + Q^ n >y + d^ n y <) ' ( 187 ) 

exist solely because of entrainment. The same procedure as in the previous two examples leads 
to the dispersion relation 

(b» [Bgo + AZ] ^r) (b s [B s 00 + A^ ] ^) - (-4 ns - * n s - ^o) 2 = • (188) 

This is quadratic in A^A^/Zcq meaning there are two sound speeds. This is a natural result of the 
doubling of fluid degrees of freedom. 

The sound speed analysis is local, but its results are seen globally in the analysis of modes of 
oscillation of a fluid body. For a neutron star, the full spectrum of modes is quite impressive 
(see McDermott et al |Z2]): polar (or spheroidal) f-, p-, and g-modes and the axial (or toroidal) 
r-modes. Epstein [HH] was the first to suggest that there should be even more modes in supernuid 
neutron stars because the superfluidity allows the neutrons to move independently of the protons. 
Mendell [73] developed this idea further by using an analogy with coupled pendulums. He argued 
that the new modes should feature a counter-motion between the neutrons and protons, i.e. as the 
neutrons move out radially, say, the protons will move in. This is in contrast to ordinary fluid 
motion that would have the neutrons and protons move in more or less "lock-step" . Analytical and 
numerical studies 1691 1351 |S] have confirmed this basic picture and the new modes of oscillation 
are known as supernuid modes. 
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12 The Newtonian limit and the Euler equations 



One reason relativistic fluids are needed is that they can be used to model neutron stars. However, 
even though neutron stars are clearly general relativistic objects, one often starts with good old 
Newtonian physics when considering new applications. The main reason is that effects (such as 
acoustic modes of oscillation) which are primarily due to the fluids that make up the star can often 
be understood qualitatively from just Newtonian calculations. There are also certain regimes in 
a neutron star where the Newtonian limit is sufficient quantitatively (such as the outer layers). 

There has been much progress recently in the analysis of Newtonian multiple fluid systems. Prix 
86 has developed an action-based formalism, analogous to what has been used here. Carter and 
Chamel ^|]|2U1|^ have done the same except that they use a fully spacetime covariant formalism. 
We will be somewhat less ambitious (for example, as in 0) by extracting the Newtonian equations 
as the non-relativistic limit of the fully relativistic equations. Given the results of Prix and Carter 
and Chamel we can think of this exercise as a consistency check of our equations. 

The Newtonian limit consists of writing the general relativistic field equations to order c° where 
c is the speed of light. The Newtonian equations are obtained in the limit that c — + oo. The 
gravitational potential, denoted <E>, is assumed to be small in the sense that 



i « 4 < ° • ( i89 ) 



To order c° the metric becomes 



ds 2 = -c 2 ^1 + ^Jj dt 2 + g ij dx i dx j , (190) 

where the x % are Cartesian coordinates and gy is the metric of Eq. (|23|) . 

If t is the proper time as measured along a fluid element worldline, then the curve it traces out 
can be written 

x^r) = {t(r),x i (r)} . (191) 
Recall that its four- velocity is given by (cf. Eq. (JSJ) 

dr^ 

u" = ^- , (192) 
dr 

but because the speed of light has been restored, u^u^ — — c 2 . The proper time is now given by 

- d,s 2 = c 2 dr 2 = c 2 (l + ^ - n£L\ dt 2 , (193) 

where v % = dx l /dt is the Newtonian three- velocity of the fluid. It is considered to be small in the 
sense that 

-L-L << 1 . (194) 
c 

Hence, to the correct order the four-velocity components are 

Mt = 1 -| + ^> M *^> ( 195 ) 
where v 2 — gijV l vK Note that the particle number current becomes 

^niu^/c) . (196) 
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BE 

A» = IT ■ ( 20 °) 



In order to write the single fluid Euler equations, keeping terms to the required order, it is 
necessary to explicitly break up the "master" function into its mass part and internal energy part 
£, i.e., to write A as 

A = -mnc 2 - £{n) , (197) 
where m is the particle mass. The internal energy is small compared to the mass, meaning 

0<— — «1. (198) 
mnc 2 

With this choice, a variation of A that leaves the metric fixed yields: 

dA = - (mc 2 + a«) dn , (199) 

where 

// = 

dn 

The B n coefficient reduces to 

B n = mc 2 + H . (201) 
n 

In terms of these variables, the pressure p is seen to be 

p = -£ + fin . (202) 

The Newtonian limit of the general relativistic fluid equations, in a general coordinate basis, is 
then 

= — + Vi(m/) , (203) 

and 

= JU* + v^VjV 1 + g^Vj + > ( 204 ) 

where Vj is the Euclidean space covariant derivative. Of course, the gravitational potential <E> is 
determined by a Poisson equation with mn as source. 

A Newtonian two-fluid system can be obtained in a similar fashion. As discussed in Sec. 
the main difference is that we obviously need two sets of worldlines, describable, say, by curves 
x£(t x ) where t x is the proper time along a constituent's worldline. Of course, entrainment also 
comes into play. Its presence implies that the relative flow of the fluids is required to specify 
the local thermodynamic state of the system, and that the momentum of a given fluid is not 
simply proportional to that fluid's flux. This is the situation for superfluid He 4 IHI3 1104) , where 
the entropy can flow independently of the superfluid Helium atoms. Superfluid He 3 can also be 
included in the mixture, in which case there will be a relative flow of the He 3 isotope with respect 
to He 4 , and relative flows of each with respect to the entropy |107) . 

Let us consider a two- fluid model like a mixture of He 4 and He 3 , or neutrons and protons in 
a neutron star. We will denote the two fluids as fluids a and b. The magnitude squared of the 
difference of three- velocities 

wV = «f - vf (205) 
is denoted w 2 so that the equation of state takes the form £ = £ (n a , rib, w 2 ). Hence, 

d£ = /i a dn a + u b dn b + adw 2 , (206) 

where 

a d£ h d£ d£ 

d™a „_,„ 2 dn h n , dw 2 n __ mv 
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The a coefficient reflects the effect of entrainment on the equation of state. Similarly, entrainmcnt 
causes the fluid momenta to be modified to 

4=< + 2-<. (208) 
The number density of each fluid obeys a continuity equation: 



-^ +Vj (n x <)=0. (209) 



Each fluid is also seen to satisfy an Euler-type equation, which ensures the conservation of total 
momentum. This equation can be written 

d + K + £x^n + V *( $ + A*) + e x icf Vit^ = , (210) 



dt 



where, 



Ax = ^ , (211) 
m 

i.e. /i x is the relevant chemical potential per unit mass, and the entrainment is included via the 
coefficients 

£x = 2p x a . (212) 
For a detailed discussion of these equations, see |8fil . 
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13 The CFS instability 



Investigations into the stability properties of rotating self-gravitating bodies are of obvious rele- 
vance to astrophysics. By improving our understanding of the relevant issues we can hope to 
shed light on the nature of the various dynamical and secular instabilities that may govern the 
spin-evolution of rotating stars. The relevance of such knowledge for neutron star astrophysics 
may be highly significant, especially since instabilities may lead to detectable gravitational- wave 
signals. In this section we will review the Lagrangian perturbation framework developed by Fried- 
man and Schutz ^2 02] for rotating non-relativistic stars. This will lead to criteria that can 
be used to decide when the oscillations of a rotating neutron star are unstable. We provide an 
explicit example proving the instability of the so-called r-modes at all rotation rates in a perfect 
fluid star. 



13.1 Lagrangian perturbation theory 

Following |41II42| . we work with Lagrangian variations. We have already seen that the Lagrangian 
perturbation AQ of a quantity Q is related to the Eulerian variation SQ by 

AQ = 6Q + C e Q (213) 

where (as before) £j is the Lie derivative that was introduced in Sec. |2J The Lagrangian change 
in the fluid velocity follows from the Newtonian limit of Eq. (|125(l : 

Av l = dtC (214) 

where £ J is the Lagrangian displacement. Given this, and 

Affy = + Vjfc (215) 

we have 

Av, = $6 + w>V<& + viVjti . (216) 

Let us consider the simplest case, namely a barotropic ordinary fluid for which £ = £{n). Then 
we want to perturb the continuity and Euler equations from the previous section. The conservation 
of mass for the perturbations follows immediately from the Newtonian limits of Eqs. i|124f> and 
(|128fl (which recall automatically satisfy the continuity equation): 

An = -nViC , tfn = -V;(nf). (217) 

Consequently, the perturbed gravitational potential follows from 

V 2 <5$ = AnGm5n = -47tGtoV j «') . (218) 

In order to perturb the Euler equations we first rewrite Eq. (|2()4(l as 

(fy+CJvi + ViU + Q-hA =0, (219) 

where ft — fi/m. This form is particularly useful since the Lagrangian variation commutes with 
the operator d t + C v - Perturbing Eq. Ij219(l we thus have 



(d t + C v )Avi + Vi [A/2 + A$ - ^A(i; 2 )) = . (220) 
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We want to rewrite this equation in terms of the displacement vector £. After some algebra 
one finds 

nfi. + 2v j V j dt^ i + (vlVjfti + V>!> + tfViVj® 

-(V i ^')V i A + V<AA = 0. (221) 

Finally, we need 

^) V 4 K 4 )+fV,A ■ (222) 
on / 

Given this, we have arrived at the following form for the perturbed Euler equation 

= . (223) 



<9n 

This equation should be compared to Eq. (15) of |4*T] . 

Having derived the perturbed Euler equations, we are interested in constructing conserved 
quantities that can be used to assess the stability of the system. To do this, we first multiply 
Eq. (|223|) by the number density n, and then write the result (schematically) as 

Ad 2 t S, + Bd t S, + C£ = , (224) 

omitting the indices since there is little risk of confusion. Defining the inner product 

(*f,&}= / V^tidV (225) 



where the asterisk denotes complex conjugation, one can now show that 

(77, AO = (£, A V )* and ( v , B£) = - <£, B v )* . (226) 

The latter requires the background relation Vj(n«*) = 0, and holds provided that n — ► at the 
surface of the star. A slightly more involved calculation leads to 

(V,C0 = {Z,C V )* . (227) 

Inspired by the fact that the momentum conjugate to is p(pt + u-' Vj)^ 1 , we now consider the 
symplectic structure 

Wfa = ^, AS t £ + ^Bf \ - (Ad t r, + ^Bn, (228) 

where r\ and £ both solve the perturbed Euler equation. Given this, it is straightforward to show 
that W{r],0) is conserved, i.e. dtW = 0. This leads us to define the canonical energy of the system 
as 

777 777 

Ec = -jW(d t £,0 = -{(d t !i,Ad t + (£,C£» . (229) 
After some manipulations, we arrive at the following explicit expression 

e c = \J { P \d t e - p\v^m 2 + pez j *ViVj(fi + $) 

\Sn\ 2 -^\^AdV (230) 
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which can be compared to Eq. (45) of |41j . In the case of an axisymmetric system, eg. a rotating 
star, we can also define a canonical angular momentum as 



Re (d v Z,A8tZ+-BZ 



(231) 



The proof that this quantity is conserved relies on the fact that (i) W(r], £) is conserved for any 
two solutions to the perturbed Euler equations, and (ii) d v commutes with pv J V j in axisymmetry, 
which means that if £ solves the Euler equations then so does <9 V £. 

As discussed by ^2 22] > the stability analysis is complicated by the presence of so-called 
"trivial" displacements. These trivials can be thought of as representing a relabeling of the 
physical fluid elements. A trivial displacement £ l leaves the physical quantities unchanged, i.e. is 
such that Sn = 8v l = 0. This means that we must have 



where, in order to satisfy the second equations, the vector \k must have time-dependence such 
that 



This means that the trivial displacement will remain constant along the background fluid trajec- 
tories. Or, as Friedman and Schutz |41j put it, the "initial relabeling is carried along with the 
unperturbed motion." 

The trivials have the potential to cause trouble because they affect the canonical energy. Before 
one can use the canonical energy to assess the stability of a rotating configuration one must 
deal with this "gauge problem". To do this one should ensure that the displacement vector £ 
is orthogonal to all trivials. A prescription for this is provided by Friedman and Schutz |41| . 
In particular, they show that the required canonical perturbations preserve the vorticity of the 
individual fluid elements. Most importantly, one can also prove that a normal mode solution is 
orthogonal to the trivials. Thus, normal mode solutions can serve as canonical initial data, and 
be used to assess stability. 

13.2 Instabilities of rotating perfect fluid stars 

The importance of the canonical energy stems from the fact that it can be used to test the stability 
of the system. In particular, 

• dynamical instabilities are only possible for motions such that E c = 0. This makes intuitive 
sense since the amplitude of a mode for which E c vanishes can grow without bounds and still 
obey the conservation laws. 

• if the system is coupled to radiation (eg. gravitational waves) which carries positive energy 
away from the system (which should be taken to mean that dtE c < 0) then any initial data 
for which E c < will lead to an unstable evolution. 

Consider a real frequency normal- mode solution to the perturbation equations, a solution of 
form £ = £ e «M+ m ¥>) _ One can readily show that the associated canonical energy becomes 




(232) 
(233) 



(234) 



(d t + C v ) X k = . 



(235) 



(236) 
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where the expression in the bracket is real. For the canonical angular momentum we get 



Jc 



(237) 



Combining Eq. (|236|l and (|237|) we see that, for real frequency modes, we will have 



E c — J c — (J p J c (238) 

m 



where a p is the pattern speed of the mode. 

Now notice that Eq. (|237(l can be rewritten as 



= — mui + m- — — a — ■ (239) 



Using cylindrical coordinates, and v 3 = flip 3 , one can show that 

- ip^VjC = pQ[m\i\ 2 + z(r x |)J . (240) 

But 

\(txi) z \<\£\ 2 (241) 

and we must have (for uniform rotation) 

<T P -n(i + -)< ¥^r < °« - « ( i - -) • (242) 

Eq. 1242|l forms a key part of the proof that rotating perfect fluid stars are generically unstable 
in the presence of radiation 02] • The argument goes as follows: Consider modes with finite 
frequency in the O — > limit. Then Eq. (|242l) implies that co-rotating modes (with a p > 0) must 
have J c > 0, while counter- rotating modes (for which a p < 0) will have J c < 0. In both cases 
E c > 0, which means that both classes of modes are stable. Now consider a small region near a 
point where a p = (at a finite rotation rate). Typically, this corresponds to a point where the 
initially counter-rotating mode becomes co-rotating. In this region J c < 0. However, E c will 
change sign at the point where a p (or, equivalently, the frequency lj) vanishes. Since the mode 
was stable in the non-rotating limit this change of sign indicates the onset of instability at a critical 
rate of rotation. 



13.3 The r-mode instability 

In order to further demonstrate the usefulness of the canonical energy, let us prove the instability 
of the inertial r-modes. For a general inertial mode we have (cf. |7U| who provide a discussion of 
the single fluid problem using notation which closely resembles the one we adopt here) 

v l ~ 8v i - f ~ H and 6$ - Sn - fl 2 . (243) 

If we also assume axial-led modes, like the r-modes, then we have 5v r ~ ft 2 and the continuity 
equation leads to 

VM - ft 3 -> ViC ~ £l 2 . (244) 
Under these assumptions we find that E c becomes (to order VI 2 ) 

E *~\j P [l^l 2 " l«*V*e| a + V 4 V, (* + £)] dV . (245) 
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We can rewrite the last term using the equation governing the axisymmetric equilibrium. Keeping 
only terms of order il 2 we have 



fT'ViVj ($ + A) ~ ItfteViVjir 2 sin 2 ( 



A bit more work then leads to 
1 



pi 21 



and 



0'^"A'V,V,(r- sin 2 0) = fiV [cos 2 6>|£ e | 2 + sin 2 6\( 



^Vi^| 2 = ft 2 {m 2 |£| 2 - 2imr 2 sin6>cos6> - f^ 9 *] 



[cos 2 0|e e 



sin 



2 «h 2 ]} 



which means that the canonical energy can be written in the form 



p {(mil - u)(mfi + 



2imil 2 r 2 sin6>cos6> - f*^ ] } 



for an axial-led mode. 

Introducing the axial stream function {/ we have 
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^ _ 
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zC7 



where Y™ — Y l m (9, <p) are the standard spherical harmonics. This leads to 



\U\< 



sm 



and 



ir 2 sin 6» cos 



1 COSfl 



m\U\ 2 [Y^dgYT* + Y^dgY^} 



After performing the angular integrals, we find that 



E r = - 



1(1+1) 



(mil — u>)(mil + u) 



2m 2 Q 2 
1(1 + 1) 



p\U\ 2 dr 



Combining this with the r-mode frequency |7(J| 



uj = mil 



1(1 + 1) 



(246) 
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(249) 

(250) 
(251) 

(252) 



(253) 



(254) 



(255) 



we see that E c < for all Z > 1 r-modes, i.e. they are all unstable. The I = m = 1 r-mode is a 
special case, leading to E c = 0. 
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13.4 The relativistic problem 



The theoretical framework for studying stellar stability in general relativity was mainly developed 
during the 1970s, with key contributions from Chandrasekhar and Friedman |23 ES] and Schutz 
|96l 197) . Their work extends the Newtonian analysis discussed above. There are basically 
two reasons why a relativistic analysis is more complicated than the Newtonian one. Firstly, 
the problem is algebraically more complex because one must solve the Einstein field equations 
in addition to the fluid equations of motion. Secondly, one must account for the fact that a 
general perturbation will generate gravitational waves. The work culminated in a series of papers 
|4()l 1411 ESI K-S9) in which the role that gravitational radiation plays in these problems was explained, 
and a foundation for subsequent research in this area was established. The main result was that 
gravitational radiation acts in the same way in the full theory as in a post-Newtonian analysis of 
the problem. If we consider a sequence of equilibrium models, then a mode becomes secularly 
unstable at the point where its frequency vanishes (in the incrtial frame). Most importantly, the 
proof does not require the completeness of the modes of the system. 
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14 Modelling dissipation 



Although the inviscid model provides a natural starting point for any investigation of the dynamics 
of a fluid system, the effects of various dissipative mechanisms are often key to the construction 
of a realistic model. Consider, for example, the case of neutron star oscillations and possible 
instabilities. While it is interesting from the conceptual point of view to establish that an instability 
(such as the gravitational-wave driven instability of the fundamental f-mode or the inertial r-mode 
discussed above) may be present in an ideal fluid, it is crucial to establish that the instability 
actually has opportunity to grow on a reasonably short timescale. To establish this, one must 
consider the most important damping mechanisms and work out whether they will suppress the 
instability or not. A recent discussion of these issues in the context of the r-mode instability can 
be found in 0]. 

From the point of view of relativistic fluid dynamics, it is clear already from the outset that we 
are facing a difficult problem. After all, the Fourier theory of heat conduction leads to instanta- 
neous propagation of thermal signals. The fact that this non-causality is built into the description 
is unattractive already in the context of the Navier-Stokes equations. After all, one would ex- 
pect heat to propagate at roughly the mean molecular speed in the system. For a relativistic 
description non-causal behaviour would be truly unacceptable. As work by Lindblom and Hiscock 
|5()| has established, there is a deep connection between causality, stability, and hyperbolicity of a 
dissipative model. One would expect an acceptable model to be hyperbolic, not allowing signals 
to propagate superluminally. 

Our aim in this section is to discuss the three main models that exist in the literature. We 
first consider the classic work of Eckart [HEI and Landau and Lifshitz |§2], which is based on a 
seemingly natural extension of the inviscid equations. However, a detailed analysis of Lindblom 
and Hiscock |51l 152) has demonstrated that these descriptions have serious flaws and must be 
considered unsuitable for practical use. However, having discussed these models it is relatively 
easy to extend them in the way proposed by Israel and Stewart [1011 1541 155) . Their description, 
the derivation of which was inspired by early work of Grad 03 an d Muller [23 and results from 
relativistic kinetic theory, provides a framework that is generally accepted as meeting the key 
criteria for a relativistic model f2D, ■ Finally, we describe Carter's more recent approach to the 
problem. This model is elegant because it makes maximal use of a variational principle argument. 
The construction is also more general than that of Israel and Stewart. In particular, it shows 
how one would account for several dynamically independent interpenetrating fluid species. This 
extension is important for, for example, the consideration of relativistic superfluid systems. 

14.1 The "standard" relativistic models 

It is natural to begin by recalling the equations that describe the evolution of a perfect fluid, see 
Sec. For a single particle species, we have the number current n M which satisfies 

V^r?/ 1 = , where tt/ 1 = nu^ . (256) 

The stress energy tensor 

T^v = P ±» +pu»u u , (257) 
where p and p represent the pressure and the energy density, respectively, satisfies 

V M T^ = . (258) 

In order to account for dissipation we need to introduce additional fields. First introduce a vector 
v 11 representing particle diffusion. That is, let 

n» — ► n?/ 1 + u* 1 , (259) 
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and assume that the diffusion satisfies the constraint u^v^ = 0. This simply means that it is 
purely spatial according to an observer moving with the particles in the inviscid limit, exactly 
what one would expect from a diffusive process. Next we introduce the heat flux g M and the 
viscous stress tensor, decomposed into a trace-part r (not to be confused with the proper time) 
and a trace-free bit t^ v , such that 

T» v — ► (p + t) J."" +pu»u u + q { »u v) + , (260) 

subject to the constraints 

u% = u%„ = = v, - = . (261) 

That is, both the heat flux and the trace- less part of the viscous stress tensor are purely spatial 
in the matter frame, and is also symmetric. So far, the description is quite general. The 
constraints have simply been imposed to ensure that the problem has the anticipated number of 
degrees of freedom. 

The next step is to deduce the permissible form for the additional fields from the second law of 
thermodynamics. The requirement that the total entropy must not decrease leads to the entropy 
flux s M having to be such that 

V M s M > . (262) 
Assuming that the entropy flux is a combination of all the available vectors, we have 

= su M + /V 1 - Az/ 1 , (263) 

where [3 and A are yet to be specified. It is easy to work out the divergence of this. Then using 
the component of Eq. (|258|l along m m , i.e. 



u^ u T^ = , (264) 

and the thermodynamic relation 

which follows from assuming the equation of state s = s(p, n), and we recall that x s = s/n, one 
can show that 



V M S" = gf(V M /3-i« , 'V^ + ^-i]V / ,.y" 



x s + A - L±£) VX - ^V M A - Jv M u" - yV, Ul/ . (266) 



We want to ensure that the right-hand side of this equation is positive definite (or indefinite) . An 
easy way to achieve this is to make the following identifications 

/3=1/T, (267) 

and 

A = —x s + E±r • (268) 
nl 

Here we note that A = g/nT where g is the Gibbs free energy density. We also identify 

v n = _ aT 2 ± nv VvX (269) 
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where the "diffusion coefficient" a > 0, and the projection is needed in order for the constraint 
u^v^ = to be satisfied. Furthermore, we can use 

r = -CVX (270) 

where £ > is the coefficient of bulk viscosity, and 

q» = -kT L^ (JfVvT + u a V a u^j (271) 

with k > the heat conductivity coefficient. To complete the description, we need to rewrite the 
final term in Eq. (|266|l . To do this it is useful to note that the gradient of the four- velocity can 
generally be written 

V ^u u = a pu + - _L M „ 6 + lo^v - ayU^ (272) 
where the acceleration is defined as 

<V = rfVpitu , (273) 



the expansion is 9 = V M u M and the shear is given by 



<rpu = \ (K V a u„+ L% V a uu) -^±^ 9. (274) 

Finally, the "vorticity" follows from 2 

wpu = \ {LI V Q u M - _L£ V a u v ) . (275) 

Since we want t^ u to be symmetric, trace-free, and purely spatial according to an observer moving 
along m m , it is useful to introduce the notation 

(V) = \ ±r±l (a p „ + A ap -\ L^ A lS ^j , (276) 

where A^ v is a general two-form. In the case of the gradient of the four- velocity, it is easy to show 
that this leads to 

(V M tt„) = <t m „ (277) 

and therefore it is natural to use 

r» v = -rja^ , (278) 
where rj > is the shear viscosity coefficient. Given these relations, we can write 

TV s M - + -+ V ^ V>1 4- T ^ T ^ V > n (279) 

By construction, the second law of thermodynamics is satisfied. 

The model we have written down is quite general. In particular, it is worth noticing that we 
did not yet specify the four- velocity u M . By doing this we can obtain from the above equations 
both the formulation due to Eckart and that of Landau and Lifshitz To arrive at the 

Eckart description, we associate u M with the flow of particles. Thus we take v* 1 = (or equivalently 
cr = 0). This prescription has the advantage of being easy to implement. The Landau and Lifshitz 

2 It is important to note the difference between the vorticity formed from the momentum and the "vorticity" in 
terms of the velocity. They differ because of the entrainment, and one can show that while the former is conserved 
along the flow, the latter is not. 
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model follows if we choose the four velocity to be a timelike eigenvector of the stress-energy tensor. 
From Eq. (|2(jl)|) it is easy to see that, by setting g M = we get 

u^ v = -pu u . (280) 

This is equivalent to setting n = 0. Unfortunately, these models, which have been used in most 
applications of relativistic dissipation to date, are not at all satisfactory. While they pass the key 
test set by the second law of thermodynamics, they fail several other requirements of a relativistic 
description. A detailed analysis of perturbations away from an equilibrium state |511 demonstrates 
serious pathologies. The dynamics of small perturbations tends to be dominated by rapidly 
growing instabilities. This suggests that these formulations are likely to be practically useless. 
From the mathematical point of view they are also not acceptable since, being non-hyperbolic, 
they do not admit a well-posed initial-value problem. 



14.2 The Israel-Stewart approach 

From the above discussion we know that the most obvious strategy for extending relativistic hy- 
drodynamics to include dissipative processes leads to unsatisfactory results. At first sight, this 
may seem a little bit puzzling because the approach we took is fairly general. Yet, the formulation 
suffers from pathologies. Most importantly, we have not managed to enforce causality. Let us 
now explain how this problem can be solved. 

Our previous analysis was based on the assumption that the entropy current s M can be described 
as a linear combination of the various fluxes in the system, the four- velocity u^, the heat-flux q M 
and the diffusion . In a series of papers, Israel and Stewart jldll 1541 155| contrasted this "first- 
order" theory with relativistic kinetic theory. Following early work by Miiller |7K) and connecting 
with Grad's 14-moment kinetic theory description 07], they concluded that a satisfactory model 
ought to be "second order" in the various fields. If we, for simplicity, work in the Eckart frame 
(cf. Lindblom and Hiscock |5U]) this means that we would use the Ansatz 

= + ^-^ (Por 2 + (3iq a q a + for^) u" + ^f- + . (281) 

This expression is arrived at by asking what the most general form of a vector constructed from 
all the various fields in the problem may be. Of course, we now have a number of new (so far 
unknown) parameters. The three coefficients /3q, 0i, and f3% have a thermodynamical origin, while 
the two coefficients and a\ represent the coupling between viscosity and heat flow. From the 
above expression, we see that in the frame moving with the effective entropy density is given 

by 

- = s - ^ (/3 r 2 + [3 iqa q a + (3 2 T al3 T al3 ) . (282) 

Since we want the entropy to be maximized in equilibrium, when the extra fields vanish, we should 
have [Po, /3i, fa] > 0. We also see that the entropy flux, 

■1? *" = - [(1 + aoT)g* + air""*,] , (283) 

is affected only by the parameters ao and a,\. 

Having made the assumption Eq. H281JI the rest of the calculation proceeds as in the previous 
case. Working out the divergence of the entropy current, and making use of the equations of 
motion, we arrive at 



T 



54 



— Q 



T (fan" 



T 



(l-7o)rTV,(f)-(l- 7l )Tr;V,(f) 



T 



(284) 



In this expression it should be noted that we have introduced (following Lindblom and Hiscock) 
two further parameters, 70 and 71. They are needed because, without additional assumptions it 
is not clear how the "mixed" quadratic term should be distributed. A natural way to fix these 
parameters is to appeal to the Onsager symmetry principle |55| . which leads to the mixed terms 
being distributed "equally" and hence 70 = 71 = 1/2. 

Denoting the comoving derivative by a dot, i.e. using w^V^r = f etcetera, we see that the 
second law for thermodynamics is satisfied if we choose 



T=~C 



2 M V T 



(285) 
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and 



~2?7 



T 



T 



V M u„ - a\V^q v - -yiTq^Vy + 73V , [^ a ]r„ a 



(287) 



where the angular brackets denote symmetrization as before. In these expression we have added yet 
another two terms, representing the coupling to vorticity. These bring further "free" parameters 
72 and 73. It is easy to see that we are allowed to add these terms since they do not affect the 
entropy production. In fact, a large number of similar terms may, in principle, be considered (see 
note added in proof in |5()|). The presence of coupling terms of the particular form that we have 
introduced is suggested by kinetic theory |55j . 

What is clear from these complicated expressions is that we now have evolution equations for 
the dissipative fields. Introducing characteristic "relaxation" times 



to = CPo 

the above equations can be written 



Kpi , and <2 = 2t7/?2 



t()f + T 

tl q v + <f 

t2Tu,v + Tnv 



<[-] , 

-kT±»" [...] , 



(288) 



(289) 
(290) 
(291) 



A detailed stability analysis by Hiscock and Lindblom 50 shows that the Israel-Stewart theory 
is causal for stable fluids. Then the characteristic velocities are subluminal and the equations form 
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a hyperbolic system. An interesting aspect of the analysis concerns the potentially stabilizing role 
of the extra parameters ((3q, ao, ...). Relevant discussions of the implications for the nuclear 
equation of state and the maximum mass of neutron stars have been provided by Olson and Hiscock 
|81l l?9"] . A more detailed mathematical stability analysis can be found in the work of Kreiss et al 



Although the Israel-Stewart model resolves the problems of the first-order descriptions for near 
equilibrium situations, difficult issues remain to be understood for nonlinear problems. This is 
highlighted in work by Hiscock and Lindblom and Olson and Hiscock They consider 

nonlinear heat conduction problems and show that the Israel-Stewart formulation becomes non- 
causal and unstable for sufficiently large deviations from equilibrium. The problem appears to 
be more severe in the Eckart frame [s3] than in the frame advocated by Landau and Lifshitz (80) . 
The fact that the formulation breaks down in nonlinear problems is not too surprising. After all, 
the basic foundation is a "Taylor expansion" in the various fields. However, it raises important 
questions. There are many physical situations where a reliable nonlinear model would be crucial, 
eg. heavy-ion collisions and supernova core collapse. This problem requires further thought. 

14.3 Carter's canonical framework 

The most recent attempt to construct a relativistic formalism for dissipative fluids is due to Carter 
|17| . His approach is less "phenomenological" than the ones we have considered so far in that 
it is based on making maximal use of variational principle arguments. The construction is also 
extremely general. On the one hand this makes it more complex. On the other hand this generality 
could prove useful in more complicated cases, eg. for investigations of multi- fluid dynamics and/or 
elastic media. Given the potential that this formalism has for future applications, it is worth 
discussing it in detail. 

The overall aim is to generalize the variational formulation described in Sec. [HI in such a way 
that viscous "stresses" are accounted for. Because the variational foundations are the same, the 
number currents play a central role (x is a constituent index as before). In addition we can 
introduce a number of viscosity tensors , which we assume to be symmetric (even though it is 
clear that such an assumption is not generally correct 0). The index X is "analogous" to the 
constituent index, and represents different viscosity contributions. It is introduced in recognition 
of the fact that it may be advantageous to consider different kinds of viscosity, eg. bulk- and shear 
viscosity, separately. As in the case of the constituent index, a repeated index X does not imply 
summation. 

The key quantity in the variational framework is the master function A. As it is a function of 
all the available fields, we will now have A(n£, t£ u , g^ v ). Then a general variation leads to 



Since the metric piece is treated in the same way as in the non-dissipative problem we will con- 
centrate on the case 6g^ v — here. This is, obviously, no real restriction since we already know 
how the variational principle couples in the Einstein equations in the general case. In the above 



As in the non-dissipative case, the variational framework suggests that the equations of motion 




(292) 
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can be written as a force-balance equation; 

V„T£ = + = (294) 

x S 

where the generalized forces work out to be 

ft = + <V[<rl$] (295) 

and 

ff=^,rr+r^(v a nf v -^V p ^ . (296) 
Finally, the stress-energy tensor becomes 

x E 

with the generalized pressure given by 

* = A-EmX-^E t ^- (298) 

x E 

For reasons that will become clear shortly, it is useful to introduce a set of "convection vectors" . 
In the case of the currents, these are taken as proportional to the fluxes. This means that we 
introduce 0£ according to 

h x 0» = < , = — > ^x = -MX • (299) 

With this definition we can introduce a projection operator 

±r=<r+Mx73x — J-&, =-L£" = . (300) 
From the definition of the force density we can now show that 

VX = (301) 

and 

KC^l =±l a f° (302) 

where the Lie-derivative along 0% is £ x = £p». We see that the component of the force which 
is parallel to the convection vector 0£ is associated with particle conservation. Meanwhile, the 
orthogonal component represents the "change in momentum" along 0% . 

In order to facilitate a similar decomposition for the viscous stresses, it is natural to choose the 
conduction vector as a unit null eigenvector associated with . That is, we take 

tt= /3£ = (303) 

together with 

= g^fc and u^0^ = -1 . (304) 
Introducing the projection associated with this conduction vector, 

±-^=g^ + rf u z — L%,f% = 0. (305) 

Once we have defined 0£, we can use it to reduce the degrees of freedom of the viscosity tensors. 
So far, we have only required them to be symmetric. However, in the standard case one would 
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expect a viscous tensor to have only six degrees of freedom. To ensure that this is the case in the 
present framework we introduce the following degeneracy condition 

u*t» v = . (306) 

That is, we require the viscous tensor t£ v to be purely spatial according to an observer moving 
along uf^. With these definitions one can show that 

/3££ E 7r=„ - (307) 

where £s = £/3£ is the Lie- derivative along and 

t^ett* - -2/3£/= . (308) 

Finally, let us suppose that we choose to work in a given frame, moving with four-velocity 
(and that the associated projection operator _L£). Then we can use the following decompositions 
for the conduction vectors 

/?£ = + ^) and ^ = + w g) . (309) 

We see that /x x = 1//3 X represents a chemical type potential for species x with respect to the chosen 
frame. At the same time, it is easy to sec that /i s = l//9s is a Lorentz contraction factor. Using 
the norm of we have 

/3^ = -/3|(l-4)=-l (310) 

where w| = v£v^. Thus 



= = ^l-uf, (311) 



clearly analogous to the standard Lorentz contraction formula. 

So far the construction is quite formal, but we are now set to make contact with the physics. 
First we note that the above results allow us to demonstrate that 

U "V„7* + £ ( M X V CT < + </*) + ^ + l -^r^C^%\ = . (312) 

Recall that similar results were central to expressing the second law of thermodynamics in the 
previous two sections. To see how things work out in the present formalism, let us single out the 
entropy fluid (with index s) by defining s M = and T = [i s . To simplify the final expressions it 
is also useful to assume that the remaining species are governed by conservation laws of the form 

V„fljJ = T x (313) 

subject to the constraint of total baryon conservation 

E r * = °- ( 314 ) 

Given this, and the fact that the divergence of the stress-energy tensor should vanish, we have 

TV,*" = -J2 M X r x - E <K - E + \^r^C^A . (315) 

Finally, we can bring the remaining two force contributions together by introducing the linear 
combinations 

ECX = and = ( 316 ) 
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constrained by 



£c x = £c s = i 

X X 



(317) 



Then denning 

# = # + ££/? (318) 

E 

we have 

TV M a" = - £ M x r x - £ </ M x - \ £ > . (319) 

x^s x s 

The three terms in this expression represent, respectively, the entropy increase due to i) chemical 
reactions, ii) conductivity, and iii) viscosity. The simplest way to ensure that the second law of 
thermodynamics is satisfied is to make each of the three terms positive definite. 

At this point, the general formalism must be completed by some suitably simple model for the 
various terms. A reasonable starting point would be to assume that each term is linear. For the 
chemical reactions this would mean that we expand each T x according to 

r x = -£c xy M y (320) 

where C xy is a positive definite (or indefinite) matrix composed of the various reaction rates. 
Similarly, for the conductivity term it is natural to consider "standard" resistivity such that 

# = -£*W ■ ( 321 ) 

y 

Finally, for the viscosity we can postulate a law of form 

= -t£""C s t$, , (322) 
where we would have, for an isotropic model, 

^ = m ±£ p ±t +~(rm Ce) ± s ^e CT (323) 

where the coefficients r/s and £s can be recognized as representing shear- and bulk viscosity, 
respectively. 

A detailed comparison between Carter's formalism and the Israel-Stewart framework has been 
carried out by Priou |84j . He concludes that the two models, which are both members of a 
large family of dissipative models, have essentially the same degree of generality and that they are 
equivalent in the limit of linear perturbations away from a thermal equilibrium state. Providing 
explicit relations between the main parameters in the two descriptions, he also emphasizes the key 
point that analogous parameters may not have the same physical interpretation. 

In developing his theoretical framework, Carter argued in favour of an "off the peg" model for 
heat conducting models |15|. This model is similar to that introduced in Sec. 1101 and was intended 
as a simple, easier to use alternative to the Israel-Stewart construction. In the particular example 
discussed by Carter, he chooses to set the entrainment between particles and entropy to zero. 
This was done in order to simplify the discussion. But, as a discussion by Olson and Hiscock [H2] 
shows, it has disastrous consequences. The resulting model violates causality in two simple model 
systems. As discussed by Priou [S3] and Carter and Khalatnikov [22], this breakdown emphasizes 
the importance of the entrainment effect in these problems. 
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14.4 Remaining issues 



We have discussed some of the models that have been constructed in order to incorporate dissipative 
effects in a relativistic fluid description. We have seen that the most obvious ways of doing this, 
the "text-book" approach of Eckart and Landau and Lifshitz fail completely. They 
do not respect causality and have serious stability problems. We have also described how this 
problem can be fixed by introducing additional dynamical fields. We discussed the formulations 
of Israel and Stewart [11)111541155) and Carter ^7] in detail. From our discussion it should be clear 
that these models are examples of an extremely large family of possible theories for dissipative 
relativistic fluids. Given this wealth of possibilities, can we hope to find the "right" model? To 
some extent, the answer to this question relies on the extra parameters one has introduced in the 
theory. Can they be constrained by observations? This question has been discussed by Geroch 
031 an d Lindblom The answer seems to be no, we should not expect to be able to use 

observations to single out a preferred theoretical description. The reason for this is that the 
different models relax to the Navier-Stokes form on very short timescales. Hence, one will likely 
only be able to constrain the standard shear- and bulk viscosity coefficients etc. Related questions 
concern the practicality of the different proposed schemes. To a certain extent, this is probably a 
matter of taste. Of course, it may well be that the additional parameters required in a particular 
model are easier to extract from microphysics arguments. This could make this description easier 
to use in practice, which would be strong motivation for preferring it. Of course, there is no 
guarantee that the same formulation will be ideal for all circumstances. Clearly, there is scope for 
a lot more research in this problem area. 
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15 Heavy ion collisions 



Relativistic fluid dynamics has regularly been used as a tool to model heavy ion collisions. The 
idea of using hydrodynamics to study the process of multiparticle production in high-energy hadron 
collisions can be traced back to work by, in particular, Landau in the early 1950s, see |11| . In the 
early days these phenomena were observed in cosmic rays. The idea to use hydrodynamics was 
resurrected as collider data became available ^3] and early simulations were carried out at Los 
Alamos [21 EJ- More recently, modelling has primarily been focussed on reproducing data from, 
for example, CERN. A useful review of this active area of research can be found in |3(J) . 

In a hydrodynamics based model, a high-energy nuclear collision is viewed in the following 
way: In the centre of mass frame two Lorentz contracted nuclei collide and, after a complex 
microscopic process, a hot dense plasma is formed. In the simplest description this matter is 
assumed to be in local thermal equilibrium. The initial thermalization phase is, of course, out 
of reach for hydrodynamics. In the model, the state of matter is simply specified by the initial 
conditions, eg. in terms of distributions of fluid velocities and thcrmodynamical quantities. Then 
follows a hydrodynamical expansion, which is described by the standard conservation equations for 
energy/momentum, baryon number and other conserved quantities, such as strangeness, isotopic 
spin, etc. (see |37| for a variational principle derivation of these equations). As the expansion 
proceeds, the fluid cools and becomes increasingly rarefied. This eventually leads to the decoupling 
of the constituent particles, which then do not interact until they reach the detector. 

Fluid dynamics provides a well defined framework for studying the stages during which matter 
becomes highly excited and compressed and, later, expands and cools down. In the final stage 
when the nuclear matter is so dilute that nucleon-nucleon collisions are infrequent, hydrodynamics 
ceases to be valid. At this point additional assumptions are necessary to predict the number of 
particles, and their energies, which may be formed (to be compared to data obtained from the 
detector). These are often referred to as the "freeze-out" conditions. The problem is complicated 
by the fact that the "freeze-out" typically occurs at a different time for each fluid cell. 

Even though the application of hydrodynamics in this area has led to useful results, it is 
clear that the theoretical foundation for this description is not a trivial matter. Basically, the 
criteria required for the equations of hydrodynamics to be valid are 1) many degrees of freedom 
in the system, 2) a short mean free path, 3) a short mean stopping length, 4) a sufficient reaction 
time for thermal equilibration, and 5) a short de Broglie wavelength (so that quantum mechanics 
can be ignored). An interesting aspect of the hydrodynamic models is that they make use of 
concepts largely outside traditional nuclear physics: e.g., thermodynamics, statistical mechanics, 
fluid dynamics, and of course elementary particle physics. This is natural since the very hot, 
highly excited matter has a large number of degrees of freedom. But it is also a reflection of the 
basic lack of knowledge. As the key dynamics is uncertain, it is comforting to resort to standard 
principles of physics like the conservation of momentum and energy. 

Another key reason why hydrodynamic models are favoured is the simplicity of the input. Apart 
from the initial conditions which specify the masses and velocities, one needs only an equation of 
state and an Ansatz for the thermal degrees of freedom. If one includes dissipation one must in 
addition specify the form and magnitude of the viscosity and heat conduction. The fundamental 
conservation laws are incorporated into the Euler equations. In return for this relatively modest 
amount of input, one obtains the differential cross sections of all the final particles, the composition 
of clusters, etc. Of course, before one can confront the experimental data, one must make additional 
assumptions about the freeze-out, chemistry, etc. A clear disadvantage of the hydrodynamics 
model is that much of the microscopic dynamics is lost. 

Let us discuss some specific aspects of the hydrodynamics that has been used in this area. As 
we will recognize, the issues that need to be addressed for heavy-ion collisions are very similar to 
those faced in studies of relativistic dissipation theory and multi-fluid modelling. The one key 
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difference is that the problem only requires special relativity, there is no need to worry about the 
spacetime geometry. Of course, it is still convenient to use a fully covariant description since one 
is then not tied down to the use of a particular set of coordinates. 

In many studies of heavy ions a particular frame of reference is chosen. As we know from our 
discussion of dissipation and causality, see Sec. 1141 this is an issue that must be approached with 
some care. In the context of heavy-ion collisions it is common to choose u M as the velocity of 
either energy transport (the Landau-Lifshitz frame) or particle transport (the Eckart frame) . It is 
recognized that the Eckart formulation is somewhat easier to use and that one can let u M be either 
the velocity of nucleon or baryon number transport. On the other hand, there are cases where the 
Landau-Lifshitz picture has been viewed as more appropriate. For instance, when ultrarelativistic 
nuclei collide they virtually pass through one another leaving the vacuum between them in a highly 
excited state causing the creation of numerous particle- antiparticle pairs. Since the net baryon 
number in this region vanishes, the Eckart definition of the four- velocity cannot easily be employed. 
This discussion is a clear reminder of the situation for viscosity in relativity, and the resolution is 
likely the same. A true frame-independent description will need to include several distinct fluid 
components. 

Multi-fluid models have, in fact, often been considered for heavy-ion collisions. One can, for 
example, treat the target and projectile nuclei as separate fluids to admit interpenetration, thus 
arriving at a two-fluid model. One could also use a relativistic multi-fluid model to allow for 
different species, e.g. nucleons, deltas, hyperons, pions, kaons, etc. Such a model could account 
for the varying dynamics of the different species and their mutual diffusion and chemical reactions. 
The derivation of such a model would follow closely our discussion in Sec. ED In the heavy- 
ion community, it has been common to confuse the issue somewhat by insisting on choosing a 
particular local rest frame at each space-time point. This is, of course, complicated since the 
different fluids move at different speeds relative to any given frame. For the purpose of studying 
heavy ion collisions in the baryon-rich regions of space, the standard option seems to be to define 
the "baryonic Lorentz frame" . This is the local Lorentz frame in which the motion of the center- 
of-baryon number (analogous to the center-of-mass) vanishes. 

The main problem with the one-fluid hydrodynamics model is the requirement of thermal 
equilibrium. In the hydrodynamic equations of motion it is implicitly assumed that local thermal 
equilibrium is "imposed" via an equation of state of the matter. This means that the relaxation 
timescale and the mean-free path should be much smaller than both the hydrodynamical timescale 
and the spatial size of the system. It seems reasonable to wonder if these conditions can be 
met for hadronic and nuclear collisions. On the other hand, from the kinematical point of view, 
apart from the use of the equation of state, the equations of hydrodynamics are nothing but 
conservation laws of energy and momentum, together with other conserved quantities such as 
charge. In this sense, for any process where the dynamics of flow is an important factor, a 
hydrodynamic framework is a natural first step. The effects of a finite relaxation time and mean- 
free path might be implemented at a later stage by using an effective equation of state, incorporating 
viscosity and heat conductivity, or some simplified transport equations. This does, of course, lead 
us back to the challenging problem of designing a causal relativistic theory for dissipation, see 
Sec. EI In the context of heavy-ion collisions no calculations have yet been performed using a 
fully three-dimensional, relativistic theory which includes dissipation. In fact, considering the 
obvious importance of entropy, it is surprising that so few calculations have been reported for 
either relativistic or nonrelativistic hydrodynamics (although see |56|). An interesting comparison 
of different dissipative formulations is provided in |77II78) . 
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16 Superfluids and broken symmetries 



In this section we discuss models that result when additional constraints are made on the prop- 
erties of the fluids. We focus on the modelling of superfluid systems, using as our example the 
case of superfluid He 4 . The equations describing more complex systems are readily obtained by 
generalizing our discussion. We contrast three different descriptions: (i) the model that follows 
from the variational framework that has been our prime focus so far if we impose that one con- 
stituent is irrotational, (ii) the potential formulation due to Khalatnikov and collaborators, and 
(iii) a "hybrid" formulation which has been used in studies of heavy-ion collisions with broken 
symmetries. 

16.1 Superfluids 

Neutron star physics provides ample motivation for the need to develop a relativistic description 
of superfluid systems. As the typical core temperatures (below 10 s K) are far below the Fermi 
temperature of the various constituents (of the order of 10 12 K for baryons) neutron stars are 
extremely cold on the nuclear temperature scale. This means that, just like ordinary matter at 
near absolute zero temperature, the matter in the star will most likely freeze to a solid or become 
superfluid. While the outer parts of the star, the so-called crust, form an elastic lattice the inner 
parts of the star are expected to be superfluid. In practice, this means that we must be able 
to model mixtures of superfluid neutrons and superconducting protons. It is also likely that we 
need to understand superfluid hyperons and colour superconducting quarks. There are many hard 
physics questions that need to be considered if we are to make progress in this area. In particular, 
we need to make contact with microphysics calculations that determine the various parameters 
of such multi-fluid systems. However, we will ignore this aspect and focus on the various fluid 
models that have been used to describe relativistic superfluids. 

One of the key features of a pure superfluid is that it is irrotational. Bulk rotation is mimicked 
by the formation of vortices, slim "tornadoes" representing regions where the superfluid degeneracy 
is broken. In practice, this means that one would often, eg. when modelling global neutron star 
oscillations, consider a macroscopic model where one "averages" over a large number of vortices. 
The resulting model would closely resemble the standard fluid model. Of course, it is important 
to remember that the vortices are present on the microscopic scale, and that they may affect 
the values of various parameters in the problem. There are also unique effects that are due to 
the vortices, eg. the mutual friction that is thought to be the key agent that counteracts relative 
rotation between the neutrons and protons in a superfluid neutron star core |74) . 

For the present discussion, let us focus on the simplest model problem of superfluid He 4 . We 
then have two fluids, the superfluid Helium atoms with particle number density n n and the entropy 
with particle number density n s . From the derivation in Sec. El we know that the equations of 
motion can be written 

V M < = (324) 

and 

ngV^/iJ] = • (325) 

To make contact with other discussions of the superfluid problem |22l 1231 1241 126) . we will use the 
notation = n£ and 0„ = [i s v . Then the equations that govern the motion of the entropy 
obviously become 

V M s^ = and s p V [p Q v] = . (326) 
Now, since the superfluid constituent is irrotational we will have 

V [( X] = . (327) 
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Hence, the second equation of motion is automatically satisfied once we impose that the fluid is 
irrotational. The particle conservation law is, of course, unaffected. This example shows how easy 
it is to specify the equations that we derived earlier to the case when one (or several) components 
are irrotational. It is worth emphasizing that it is the momentum that is quantized, not the 
velocity. 

It is instructive to contrast this description with the potential formulation due to Khalatnikov 
and colleagues |58l 164) . We can obtain this alternative formulation in the following way |23| . 
First of all, we know that the irrotationality condition implies that the particle momentum can be 
written as a gradient of a scalar potential a (say). That is, we have 

v p = = -V p a . (328) 
m 

Here m is the mass of the Helium atom and v p is what is traditionally (and somewhat confusedly) 
referred to as the "superfluid velocity". We see that it is really a rescaled momentum. Next 
assume that the momentum of the remaining fluid (in this case, the entropy) is written 

H* = <d P = k p + V p . (329) 

Here n p is Lie transported along the flow provided that s p n p — (assuming that the equation of 
motion Eq. I)326[) is satisfied). This leads to 

s p V p (/> = s p e p . (330) 

There is now no loss of generality in introducing further scalar potentials (3 and 7 such that 
K p = /3V P 7, where the potentials are constant along the flow-lines as long as 

S PV P P = s p V P 7 = . (331) 

Given this, we have 

6 P = V P 4> + PV p1 . (332) 

Finally, comparing to Khalatnikov's formulation [SSJ |M] we define O p = —kw p and let cj> — > nC, 
and (3 — > n[3. Then we arrive at the final equation of motion 

- — = w p = -V P C - /3V p7 . (333) 

Eqs. 1)328)1 and 1)333)1 . together with the standard particle conservation laws, are the key equations 
in the potential formulation. As we have seen, the content of this description is identical to that 
of the canonical variational picture that we have focussed on in this review. We have also seen 
how the various quantities can be related. Of course, one has to exercise some care in using this 
description. After all, referring to the rescaled momentum as the "superfluid velocity" is clearly 
misleading when the entrainment effect is in action. 



16.2 Broken symmetries 

In the context of heavy-ion collisions, models accounting for broken symmetries have sometimes 
been considered. At a very basic level, a model with a broken U(l) symmetry should correspond 
to the superfluid model described above. However, at first sight our equations differ from those 
used in, for example, |l()()l 1871 1118) . Since we are keen to convince the reader that the variational 
framework we have discussed in this article is able to cover all cases of interest (in fact, we believe 
that it is more powerful than alternative formulations) a demonstration that we can reformulate 
our equations to get those written down for a system with a broken U(l) symmetry has some merit. 
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The exercise is also of interest since it connects with models that have been used to describe other 
superfluid systems. 

Take as the starting point the general two-fluid system. From the discussion in Section ITU1 we 
know that the momenta are in general related to the fluxes via 

H* = B x n* + A xy n y p . (334) 

Suppose that, instead of using the fluxes as our key variables, we want to consider a "hybrid" 
formulation based on a mixture of fluxes and momenta. In the case of the particle-entropy system 
discussed in the previous section, we can then use 

1 A ns 

< = w y; - ■ (335) 

Let us impose irrotationality on the fluid by representing the momentum as the gradient of a scalar 
potential <p. With p? = V p <^ we get 

1 A na 

< = rrVpf ~ ■ (336) 

B n B n 

Now take the preferred frame to be that associated with the entropy flow, i.e. introduce the unit 
four velocity u p such that = n s u p = su p . Then we have 

n n p = nu p - V 2 V p ip (337) 

where we have defined 

n = -— — and V 2 = -— . (338) 

B n B n v ' 

With these definitions, the particle conservation law becomes 

V p < = V p (nu p - V 2 V p Lp) = . (339) 

The chemical potential in the entropy frame follows from 

jj, = -u p n n p = -u p V p ip . (340) 

One can also show that the stress-energy tensor becomes 

T p v = + (* + p)u p u v - V 2 W p (fW^ (341) 

where the generalized pressure is given by as usual, and we have introduced 

* + p = B s s 2 + A sn sn . (342) 

The equations of motion can now be obtained from V M T M ^ = 0. (Keep in mind that the equation 
of motion for x = n is automatically satisfied once we impose irrotationality, as in the previous 
section.) This essentially completes the set of equations written down by, for example, Son [100] . 
The argument in favour of this formulation is that it is close to the microphysics calculations, which 
means that the parameters may be relatively straightforward to obtain. Against the description 
is the fact that it is a, not very elegant, hybrid where the inherent symmetry amongst the different 
constituents is lost, and there is also a risk of confusion since one is treating a momentum as if it 
were a velocity. 
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17 Final remarks 



In writing this review, we have tried to discuss the different building blocks that arc needed if one 
wants to construct a rclativistic theory for fluids. Although there are numerous alternatives, we 
opted to base our discussion of the fluid equations of motion on the variational approach pioneered 
by Taub |102j and in recent years developed considerably by Carter ^] ^5] ^]. This is an 
appealing strategy because it leads to a natural formulation for multi- fluid problems. Having 
developed the variational framework, we discussed applications. Here we had to decide what to 
include and what to leave out. Our decisions were not based on any particular logic, we simply 
included topics that were either familiar to us, or interested us at the time. That may seem a little 
peculiar, but one should keep in mind that this is a "living" review. Our intention is to add further 
applications when the article is updated. On the formal side, we could consider how one accounts 
for elastic media and magnetic fields, as well as technical issues concerning relativistic vortices 
(and cosmic strings). On the application side, we may discuss many issues for astrophysical fluid 
flows (like supernova core collapse, jets, gamma-ray bursts, and cosmology). 

In updating this review we will obviously also correct the mistakes that are sure to be found by 
helpful colleagues. We look forward to receiving any comments on this review. After all, fluids 
describe physics at many different scales and we clearly have a lot of physics to learn. The only 
thing that is certain is that we will enjoy the learning process! 
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